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PREFACE 


This  report  describes  the  finite  element  computer  program  SOE.STRUCT  used  in 
the  evaluation  of  soil-structure  interaction  of  earth  retaining  structures.  The  initial  versicm  of 
tire  program  was  developed  by  Professors  G.  W.  Clou^  and  J.  M.  Duncan  in  1969  and  has 
been  enhanced  during  the  last  20  years  by  Professor  Qougb  and  his  cowoikers.  This  report 
documents  the  version  of  the  program  in  which  a  procedure  for  performing  soil-structure 
interaction  analysis  using  die  coqded-boundary  elernem/finite  element  method  is  incorporated. 
The  work  was  performed  by  Mark  Zaico.  Ph.D..  student  at  Virginia  Tedi.  unda  the  guidance 
of  Professor  T.  Ktqrpusamy.  Professor  of  Qvil  Engineering.  Virginia  Tedi.  This  work  was 
funded  by  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  under  Oontiact  No. 
DACrW39-91-C-(X)8S.  Ftinding  for  the  adiqptatitxi  of  the  program  and  documentation  was 
provided  by  Headquarters.  U.S.  Army  Corps  of  Engineers  (HQUSACE).  under  the  Computer- 
Aided  Structural  Engineering  (CASE)  Project 

This  report  was  prepared  Dr.  T.  Kuppusamy.  Dr.  Marie  A.  Zarco.  and  Dr.  Robert 
Ebeling.  Scientific  and  Rngirreering  ^rpUcadoos  Center.  Computer-Aided  Engiireeiing  Division 
(CAED),  Information  Technology  Ltiboratory  (TTL).  WES.  The  worit  was  managed  and 
coordinated  by  Dr.  Reed  L.  Mosher.  Interdesciplinary  Research  Group.  CAED.  ITL.  and  Dr. 
Ebeling.  AU  the  work  was  acconqdidied  under  the  general  supervision  of  Dr.  Reed  Mosher. 
Acting  Chief.  CAED.  and  Dr.  N.  Radhakririman.  Director.  ITL. 

At  the  time  of  publication  of  this  report  Director  of  WES  was  Dr.  Robert  W.  Whalirt 
Commander  was  COL  Bruce  K.  Howard.  EN. 
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CONVERSION  FACTORS,  NON^  TO  SI  (METIUQ  UrarS  OF  MEASUREMENT 

Non-Sl  units  of  measurement  used  in  this  nport  can  be  converted  to  SI  (metric) 
units  as  follows: 


_ Multiply 

cubic  feet 
foet 

pounds  (mass)  per 
cubic  foot 

pounds  (force)  per 
square  foot 

pounds  (force)  per 
square  inch 

square  feet 

square  inches 


_ 

0.283 1685 
0.3048 
16.01846 

47.88026 

6.894757 

0.09290304 

6.4516 


To  Obtain 
cubic  metres 
metres 

kilograms  per 
cuHc  metre 

pascals 

Idlopascak 

square  metres 
square  centimeters 
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USER'S  GUIDE  FOR  THE  INCREMENTAL 
CONSTRUCTION  SOIL-STRUCTURE 
INTERACTION  PROGRAM  SOILSTRUCT 

WITH 

FAR-FIELD  BOUNDARY  ELEMENTS 


PART  I:  DESCRIPnON  OF  PROGRAM  SOILSTRUCT 


latrodactiOB 

1.  SOILSTRUCT  is  a  general  purpose  finite  dement  jMrogram  for  two-dimensional 
plane  strain  analysis  of  sod-structure  interaction  problems.  It  calculates  diq)lacemaits 
and  stresses  due  to  incremental  construction  andAn'  load  aj^dication  and  is  oqMble  of 
modding  nonlinear  stress-strain  material  bdnvior.  The  simulation  of  inaemental 
consmiction  may  include  end)ankment  constructkm  or  backfilling,  excavation, 
installation  of  a  ttrut  or  tie-back  anchor  excavation  support  system,  removal  of  the 
same  system,  and  the  placement  of  concrete  or  other  constructitm  materials.  The 
incronental  loading  simulation  m^  consist  of  the  applicatkm  of  concentrated  loads, 
boundary  pressures,  or  loads  due  to  ten^Miiture  chMoges  in  non-soU  materials. 

2.  The  initial  version  of  SOILSTRUCT  was  developed  by  Professors  G.  W.  Clough 
and  J.  M.  Duncan  fiar  use  in  the  analyds  of  Port  Alto  and  Old  River  U-firante  lodes 
(Clough  and  Duncan  1969).  Since  th^  the  program  has  been  modified  to  eaqpand  the 
crqiabilities  of  the  finite  demem  constitutive  models,  load  vector  formulation  algo¬ 
rithm,  the  size  of  the  problem  which  nuiy  be  analyzed,  and  the  transfia’  of  input,  out¬ 
put,  restart,  and  plot  data  files  by  means  of  disc  storage.  These  modifications  have 
bemi  nuule  in  conjunction  with  a  number  of  {^ejects  at  the  US  Aimy  En^neer 
Waterwi^  Experimem  Station  (Ebding  1990,  Ebdiiig,  Duncan  arxl  Gou^  1990, 
Ebding,  Clough,  Duncan  and  Brandon  1992,  R^alado,  Duncan  and  Clou^  1992, 
Ebding,  Mosher,  Abraham  and  Peters  1993).  This  version  of  the  program  reflects 
modifications  nuule  by  Mark  Zarco  under  the  supervition  of  Profisssm  T.  Kui^xisamy, 
to  incorporate  the  capability  of  modeling  fiu-fidd  effects  by  boundary  dements. 

3.  SOILSTRUCT  has  been  coded  in  FORTRAN  77  language  and  condsts  of  a 
tiuin  program  and  33  subroutines  tunned:  DETNA  INTTAL,  STRSIT,  QUAD, 
BAREL,  EXCAV,  EQNDFO,  SURFLD,  JTSTF,  SUBSIP,  JSTRES,  SEEP, 
MODCAL,  BUILD,  OPTSOL,  AUXOUT,  STRESS,  PRNCIP,  PRNTFD,  GETFJL, 
NOTENS,  UNBALS,  PRINTFO,  PRNTJDM  INTERF,  PRINTFD,  BEMSTF, 
MATRX,  FUNC,  MELAN,  SETCON,  NMATRX,  and  STFSYM.  A  user's  guide  for 
the  SOILSTRUCT  Program  can  be  fisund  in  Appoxlix  A.  Appendix  B  contains  the 
sequoice  of  operations  for  the  SOILSTRUCT  Program. 


4 


Previous  Contributioiis  to  SOILSTRUCT 


4.  Numerous  versions  of  the  SOILSTRUCT  code  have  beai  wntten  since  it  was 
introduced  by  Clough  and  Duncan  in  1969.  Most  modifications  were  used  to  solve  a 
specific  soil-structure  interaction  problem.  The  resulting  computer  code  continued  to 
be  referred  to  as  SOILSTRUCT.  Because  of  the  lack  of  documentation  within  the 
computer  code,  the  authors  of  this  report  were  not  able  to  identify  the  early 
contributions  contained  within  the  version  that  was  used  as  the  basis  for  the  computer 
code  described  in  this  report.  Since  1985,  four  aiudytical  enhancements  have  been 
made  to  the  version  of  SOILSTRUCT  described  in  this  report,  two  are  known  as  "the 
alpha  method,"  the  third  is  a  shear  stress-deformation  model  for  interfiu:e  elements, 
and  the  fourth  is  the  modeling  of  the  fiu’-field  dotruin  using  boundary  elements. 

5.  The  alpha  method  was  first  developed  for  analyzing  walls  that  are  laterally 
loaded  so  heavily  that  gaps  develop  along  the  interface  between  the  base  and  the 
foundation  (Ebeling,  Duncan  and  Clough  1990,  and  Ebeling,  Clough,  Duncan  and 
Brandon  1S^2).  This  numerical  method  was  later  extended  to  two-dimensional 
elements  and  used  to  reduce  numerical  inaccuracies  such  as  overshoot  stresses  in 
nonlinear  soil  elements  which  are  in  failure  or  near  failure  (Regalado,  Duncan  and 
Clough  1992).  This  situation  occurs  r^en  walls  retaining  backfills  that  are  at,  or  near, 
an  active  state  of  stress  are  analyzed. 

6.  A  third  enhancement  contained  within  this  version  of  SOILSTRUCT  is  the 
reintroduction  of  the  hyperbolic  shear-stress  displacement  relationship  for  the  interfiu:e 
element  (Ebeling,  Peters,  and  Clough  1990).  Although  present  in  the  original  version 
of  SOILSTRUCT,  this  nonlinear  displacement  relationship  was  missing  fi’om  the 
version  of  computer  code  in  vdiich  the  alpha  method  was  incorporated. 

7.  Coupling  of  the  boundary  element  method  of  analysis  with  the  finite  element 
method  is  described  in  Part  IV  of  this  report. 

Finite  Elements  and  Boundary  Elements  Employed 


8.  Three  types  of  finite  elements  are  used  to  represent  the  behavior  of  different 
nuterials:  (a)  a  two-dimen»onal  continua  demoit,  (b)  an  imeiface  element,  and  (c)  a 
one-dimenaonal  bar  element. 

9.  A  two-dimennonal,  subparametric,  quadrilateral  element  (QMS)  is  used  to  rq>re- 
sent  the  soil  and  most  structural  materials.  Structural  supports,  such  as  the  struts  or 
tidrack  components  of  an  excavation  support  ^em,  are  typically  modeled  as  a  spring 
support  using  bar  dements.  However,  two-dimendonal  dements  have  been  us^  to 
modd  these  supports.  The  geometry  of  this  element,  devdoped  by  Doherty,  \^son 
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•nd  Taylor  (1969)  is  defined  by  finir  exxemal  nodes,  v^iile  the  di^lacement  fiinctions 
include  an  intonal  fifth  node.  To  improve  flexural  re^mnse,  a  shear  itfr><in^ 

calculated  at  the  location  of  the  internal  fifth  node,  is  imposed  throughout  the 
The  QMS  element  can  be  allowed  to  degrade  to  a  triangular  element  by  lettiitg  two 
adjacent  nodes  of  the  quadrilateral  coincide. 

10.  The  Goodman,  Taylor  and  Breeke  (1968)  intcrfece  element  is  used  to  allow  for 
relative  movement  between  different  materials  such  as  between  a  soil  backfill  and  a 
support  wall.  This  element  is  defined  by  four  nodes,  with  each  of  the  two  pairs  of 
nodes  having  the  same  coordinates;  thus,  this  type  of  element  has  no  thickness. 

11.  One  dimensional,  two  node  bar  or  spring  elonents  are  used  to  model  the  behav¬ 
ior  of  a  variety  of  structural  systems.  This  includes  the  modeling  of  structural  supports 

such  as  braces  or  ti^acks  or  the  modeling  of  reinforcement  placed  within  a  soil 
backfill. 

12.  Linear  boundary  elements  based  on  the  Melan  fimdamental  solution  (Telles  and 
Br^bia  1981),  are  used  to  model  the  far-field  boundary. 


6 


PART  n:  MATERIAL  STRESS-STRAIN  BEHAVIOR 


Introductioa 

13.  Several  inodes  of  stress-strain  behavior  are  utilized  to  rqiresent  the  re^nse  of 
soil,  construction  materials,  and  the  inter&ce  re^on  between  different  materials. 

14.  The  constitutive  rdationship  used  for  all  two-dimenaonal  elements  is  Hooke's 
law.  SOILSTRUCT  uses  an  incremental,  equivalent  linear  method  of  analysis  to  model 
nonlinear  material  behavior.  In  this  type  of  analysis,  the  incremental  changes  in  stresses 
are  related  to  the  incremental  strains  through  a  linear  relationship.  This  relationship  is 
defined  for  each  structural  element  by  two  engineering  constants,  the  Young's  moduli 
and  the  Poison's  ratio. 


Nonlinear  Stress-Strain  Response  of  SoU 

15.  A  plane  strain,  isotropic  drained  or  undrained  stress-strain  soil  model  is  incor¬ 
porated  within  soil  SOILSIRUCT.  The  prognun  uses  a  nonlinear,  stress-dependent 
hyperbolic  curve  to  represent  the  relationship  between  stress-strain  response  during 
primary  loading  of  the  soil  (Figure  la)  and  a  linear  stress-strain  response  during  un¬ 
loading  or  reloading  of  the  soil  (Figure  lb).  Tlw  unload-reload  stress-strain  response  is 
applicable  when  the  current  stress  state  is  less  than  that  which  has  been  appli^  previ¬ 
ously;  otherwise,  the  primary  loading  stress-strain  is  iq>pTopriate.  Laboratory  testing 
and  interpretation  procedures  for  determining  the  parameters  used  to  define  the  soil 
model  are  described  in  Duncan,  Byrne,  Wong  and  Mabry  (1978).  A  brief  review  of  the 
hyperbolic  model  is  given  in  paragraphs  18  to  22. 

16.  The  nonlinear  soil  response  to  loading  is  modeled  by  performing  a  series  of 
analyses  in  which  each  load  is  ^iplied  incrementally,  with  the  total  change  in  stress 
computed  at  the  center  of  each  soU  element  being  equal  to  the  sum  of  the  incremental 
changes  in  stress  over  all  the  load  steps.  In  gmeral,  the  greater  the  curvature  of  the 
stress-strain  relationship  or  the  larger  the  magnitude  of  the  applied  load,  the  greater 
the  number  of  load  steps  required  to  accurately  model  the  nonlinear  soil  response.  This 
may  be  achieved  in  two  ways  using  SOILSTRUCT;  dther  the  total  load  ^proach  us¬ 
ing  a  greater  number  of  incremental  loadings,  or  during  the  course  of  each  load  case 
analysis,  the  load  vector  may  be  applied  in  a  series  of  increments  using  the  substep  op¬ 
tion. 

17.  Application  of  each  loading  in  the  finite  element  analysis  results  in  a  change  in 
stress  within  each  of  the  soil  elements.  In  addition  to  the  change  in  stress,  there  is  a 
corresponding  change  in  stififiiess.  Since  each  incremental  analysis  is  performed  assum¬ 
ing  equivalent  linear  element  response,  SOILSTRUCT  updates  the  >^ue  of  the  dastic 
moduli  assigned  to  each  soil  element.  To  account  for  the  change  in  stiffiiess  that  oc¬ 
curs  during  the  application  of  a  load  increment,  each  incremental  load  calculation  may 
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b.  Umar  vnloat£r^-rehadmg  stress-strain  relation^p 

Figure  1.  H)^>erbotic  model  for  stress-strain  behavior 
(cfkr  Duncan,  Byrne,  Wong,  cndMcUny  1978) 
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br  repeated  using  the  iteration  option.  When  the  iteration  option  is  invoked,  the  load 
vector  is  reapplied  with  a  revised  value  for  the  dement  sti^ess.  The  value  assigned 
for  the  stiffiiess  of  the  soil  element  reflect  the  average  of  the  stress  state  developing  at 
the  end  of  the  previous  load  case,  or  substep,  and  that  which  develops  during  the 
current  iteration.  However,  when  only  one  iteration  is  specified,  the  modulus  values 
are  calculated  using  the  stresses  developing  at  the  end  of  the  previous  load  increment. 
Upon  completion  of  the  last  iteration  for  each  load  case  or  substep,  the  arrays 
tabulating  the  values  of  the  total  nodal  point  displacements  and  total  element  stresses 
are  updated  with  the  computed  incremental  values. 


Primary  Loading  -  Young's  Moduli 

18.  Prior  to  each  analysis  a  tangent  Young's  modulus  is  assigned  to  each  soil 
element.  The  stress-dependent  value  of  is  computed  using  the  relationship 


(1) 


where 

Ej  =  initial  Young's  modulus 
Ry=  failure  ratio 
SL  =  stress  level 


The  initial  Young's  modulus  £/  is  equal  to 


where 

K  =  modulus  number 
Pq  =  atmospheric  pressure 
n  =  modulus  exponent 

^3  =  minor  principal  stress 


(2) 


19.  The  proportion  of  mobilized  shear  strength  for  each  soil  element  is  reflected  in 
the  value  of  the  stress  level  SL.  SL  is  equal  to  the  current  deviator  stress  (ay  -  oj) 
divided  by  the  deviator  stress  at  fiulure  (ay  -  aj)^,  denoted  by  the  subscript /. 
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(3) 


where  aj  =  major  principal  stress.  The  value  of  SL  ranges  from  a  value  equal  to  zero 
to  a  value  equal  to  unity.  SL  equal  to  zero  indicates  an  isotropic  str^  state,  >^e  SL 
equal  to  unity  corresponds  to  the  complete  mobilization  of  shear  reastance  within  the 
soil  element. 

20.  This  version  of  SOILSTRUCT  defines  the  deviator  stress  at  fiulure  using  the 
original  Duncan  formulation  (Duncan  and  Chang,  1970).  In  the  original  Duncan 
formulation,  the  value  of  the  minor  principal  stress  at  fiulure  is  set  equal  to  the  current 
minor  principal  stress.  The  de^aator  stress  at  fiulure  is  ^ven  by 


2ccos^+2a3sin^ 

l-sin<> 


(4) 


where 

c  =  cohesion  intercept 
^  =  angle  of  internal  fiiction 

21 .  The  feilure  ratio  J?^relates  the  ultimate  deviator  stress 

(<r,  -  a,)^  =  RfitTi  -  cTj)^  (5) 

The  ultimate  deviator  stress  is  the  asymptote  to  the  stress-strain  h5q)erbola,  as  in 
Figure  la.  The  value  ofi?^is  always  less  than  unity  and  varies  firom  0.5  to  0.9  for  most 
soils. 


Unload-Rdoad  Stress-Strain  Behavior  -  Young's  Modulus 

22.  During  unloading  or  reloading,  when  the  current  deviator  stress  is  less  than  that 
which  has  been  applied  during  previous  loadings,  a  stress-dependent,  linear  rwponse  is 
as  shown  in  Figure  lb.  In  this  case,  the  value  is  computed  u^g 


udiere  .ATu^  =  unload-reload  modulus  number. 
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Poison's  Rntio 


23.  The  second  elastic  parameter  used  to  d^ne  tlw  nuterial  bdiavior  of  soil  is  the 
Poison's  ratio  v.  When  using  Poison's  ratio,  two  values  are  specified;  a  constant  value 
which  is  applicable  for  all  states  of  stress  prior  to  fiulure,  5Z,  <  1,  and  the  value  of 
Poison's  ratio  applicable  \«dien  the  shear  strength  of  the  soil  is  fully  mobilized,  SL-\. 


24.  When  computing  initial  stress  by  gravity  tum-on  analysis,  the  value  of  Poison's 
ratio  used  for  the  soil  model  may  not  be  suhi^le  for  initial  stress  conditions.  For  ex¬ 
ample,  the  Poison's  ratio  for  undrained  analysis  is  generally  taken  to  be  neatly  O.S  be¬ 
cause  saturated  soil  is  near  incompressible.  By  contrast,  the  initial  stress  conditions 
should  be  fully  drained.  Therefore,  the  initial  stress  computations  are  based  on  a  value 
of  Poison's  ratio  that  gives  the  correct  ratio  for  level  ground  conditions;  that  is 

o'x/a'h  ~  ^o-  Using  the  value  of  Kq  input,  the  Poison's  ratio  used  for  initial  conditions 
is  given  by 


v= 


(7) 


Jaky  (1948)  suggested  that  Kq  may  be  i^proximated  using  the  relationship 


=  l-stn^ 

where  is  the  angle  of  internal  fiiction. 


(8) 


Ovenhoot  in  Soil  Elements  -  ALPHA  Method 

25.  The  alpha  method  is  used  to  analyze  walls  retaining  backfills  containing  soil 
elements  which  are  in  or  near  fiulure  (Reg^ado,  Duncan  and  Clough  1992).  The  alpha 
method  is  used  to  reduce  numerical  inaccurades  such  as  overshoot  stresses  in  these 
nonlinear  soil  elements. 

26.  The  finite  element  program  SOILSTRUCT  employs  hyperbolic  relationships  to 
model  the  nonlinear  stress-strain  bdiavior  of  soils.  As  discussed  predously,  nonlinear 
behavior  is  approxinuted  by  performing  the  analyses  in  a  series  of  loading  increments 
or  "steps".  >^ithin  each  step  of  the  analysis  the  materials  are  treated  as  if  their  bdmvior 
is  linear.  The  modulus  values  assigned  to  each  element  are  adjusted  in  accordance  vrith 
their  stresses  to  simulate  nonlinear  behavior. 

27.  The  way  the  stresses  in  an  element  might  change  during  one  step  of  an  analysis 
is  shown  in  Figure  2.  Figure  2(a)  shows  the  state  of  stress  at  the  center  of  an  element 
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and  the  coiTeqx)nding  Mohr's  circle  at  the  beginning  of  a  load  step.  Figure  2(b)  shows 
the  state  of  stress  and  the  Mohr's  circle  for  the  same  demoit  after  a  load  step  that 
brought  about  the  changes  in  stresses  Ac^  and  Atj^  Figure  2(b)  illustrates  that 
after  the  addition  of  the  stress  incremoits,  the  resulting  combination  of  stresses  may 
result  in  an  unacceptable  state  of  stress.  A  Mohr's  drcle  of  stress  above  the  fiulure 
envelope  is  indicative  of  overshoot. 

28.  The  ideal  situation  would  be  that  the  load  iq>plied  during  a  load  step  was  just 
large  enough  to  bring  the  most  severely  stresses  dement  to  fiulure.  This  situation 
would  correspond  to  a  Mohr's  circle  which  was  tangent  to  the  fiulure  envelope.  The 
purpose  of  the  alpha  method  is  to  achieve  reailts  closer  to  the  ideal  situation,  in  which 
no  element  is  overstressed.  The  bade  idea  behind  the  method  is  to  reduce  the 
magnitudes  of  the  loads  applied  in  any  step  of  the  analysis  to  a  point  that  no  soil 
element  incurs  a  stress  level  greater  than  unity. 

29.  Figure  3(a)  illustrates  the  same  situation  shown  in  Figure  2(b)  in  which  the 
stress  increments  during  a  load  step  resiih  in  overshoot  in  a  soil  element.  Figure  3(b) 
shows  a  Mohr's  circle  which  has  been  reduced  in  size  so  that  it  is  tangent  to  the  fiulure 
envelope.  This  reduction  can  be  accomplished  by  multiplying  the  magnitudes  of  all  the 
stress  change  components  (Ao^,  Aoy  and  Atj^)  by  a  fiictor  a  (hence  the  name  "alpha" 
method).  Reducing  the  stresses  by  this  fiictor  in  all  the  elements  in  the  mesh  can  be 
accomplished  by  reducing  the  load  by  the  same  fiictor,  since  the  system  bd".* '  i 
lineariy  during  each  step  of  the  analysis. 

30.  A  situation  may  arise  in  which  several  elements  fiul,  with  each  one  requiring  a 
dififerent  value  of  the  reduction  fiu:tor  alpha.  When  this  situation  occurs,  the  most 
severely  stressed  element  will  require  the  smallest  value  of  alpha.  Using  the  alpha  value 
to  reduce  the  load  will  ensure  that  fiulure  takes  place  only  in  the  worst  stresses  elonent 
and  nowhere  else  in  the  mesh.  This  worst  stressed  element  is  then  assigned  properties 
that  reduce  its  ability  to  carry  more  load.  After  this  change  in  the  properties  of  the 
most  severely  stressed  element,  the  remaining  portion  of  the  load  for  this  particular 
load  step  is  applied.  This  process  of  reducing  the  load  in  order  to  allow  only  one 
element  to  fiul  at  a  time  is  repeated  until  the  entire  load  for  the  load  step  has  been 
applied. 

31.  In  the  process  of  assigning  properties  to  a  fiuled  element,  consideration  is  ipvai 
to  avoiding  too  abrupt  changes  in  modulus  value.  When  an  element  first  incurs  a  stress 
level  that  indicates  ^ure,  its  modulus  is  reduced  to  one-tenth  its  value  brfore  fiulure. 
This  proc^  of  reducing  the  modulus  is  repeated  in  subsequent  loadings  for  as  long  as 
the  elemoit  is  in  a  state  of  fiulure.  In  order  for  the  modulus  value  not  to  become 
urvealistically  snudl,  the  program  user  can  specify  a  Iowa*  limit  for  the  modulus  to  be 
assigned  to  fiuled  elements.  When  the  process  of  continuous  reduction  in  the  value  of 
modulus  for  a  failed  element  results  in  a  value  smaller  than  the  lower  limit,  its  value  is 
kept  at  the  lower  limit.  It  is  also  at  this  stage  that  the  element  is  assigned  a  fiulure 


(a)  Before  Load  Step 


T 


Figure  2.  Stress  Changes  in  a  Stdl  Element  during  a  Load  Step  and  the 
Occurrence  of  Overshoot  (Regalado,  Duncan,  andChu^  1992). 
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(b)  F(Hnared  Load  end  Corre^onda^  Stress  Conqionent  Increases 

in  Soil  Element 

Figures.  Factoring  die  Load  to  Eliminate  Overshoot 
(Regalado,  Duncan  and  Chugk  1992) 
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value  of  Poisson's  fatio.  The  element  retains  these  property  values  for  as  loqg  as  it 
renuuns  in  failure. 

32.  Although  this  procedure  allows  a  gradual  transition  to  fiulure,  it  does  not 
completely  remove  the  load  carrying  capacity  of  a  fidled  elemoit.  Consequoitly,  fiuled 
elements  are  still  able  to  carry  some  additional  load,  making  it  possible  for  than  to 
exceed  a  stress  level  of  unity  by  some  usually  small  amount. 

33.  In  sumnuuy,  the  steps  in  the  alpha  method  are  as  follows:  a.  if  overshoot  occurs 
in  any  soil  element,  reduce  the  stresses  to  eliminate  ovendioot  in  the  most  severely 
stressed  element;  b.  assign  post-fiulure  properties  to  this  element;  and  c.  apply  the 
renuuning  portion  of  the  load  and  continue  the  analysis.  The  data  input  for  the  alpha 
method  are  described  in  Section  4.  of  Appendix  A.  For  additional  dettuls  regarding  the 
use  of  the  alpha  method  for  controlling  ovoshoot  in  two-dimensional  elements,  see 
Chapter  2  of  Regalado,  Duncan  and  Clough,  1992. 
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PART  m:  MODELING  STRUCTURAL  ELEMENTS 
latrodactioB 

34.  TIhs  chiqpter  desoibes  how  the  stress-strain  behavior  of  diffsieiit  structural  de¬ 
nies  used  are  modded.  Included  is  a  discusdon  on  nutterial  re^nse  of  structural 
dements  modded  usiiig  two-dimensioiud  QMS,  bar  and  interfiice  dements.  Also  dis¬ 
cussed  is  the  dpha  method  for  modetir^  base  squuation. 

Stnictnral  Material  Response 

35.  Structural  materials  such  as  wood,  concrete,  or  sted  are  modded  using  two-di- 
mensiond  QMS  dements  with  linear  elastic  stress-strain  behavior  assumed.  Support 
dements  such  as  struts  or  anchors  are  typically  modded  using  one-dimensimud  bar 
nmnbers,  and  are  also  assumed  to  behave  linea^.  Bar  dements  as  formulated  within 
in  SOILSTRUCT  have  the  capability  to  re^nd  in  compression  only,  in  tension  only, 
or  both  tension  and  conq)ression.  In  addition,  sladc  in  the  suf^xnt  system  at  the  time 
of  installation  may  be  accounted  for  by  qjedfying  an  inhid  vdue  of  di^lacement  for 
the  bar  dement. 


Interface  Response 

36.  Interfoce  dements  are  used  to  allow  for  relative  movement  between  different 
materid  regions,  such  as  between  a  soil  backfill  and  a  support  wall.  These  dement  are 
defined  by  fimr  nodes,  each  node  having  two  d^rees  of  fieedom;  each  of  the  two 
pairs  of  nodes  sharing  the  same  coordinates.  The  interfiice  dement,  therefore,  is  of 
finite  length  but  zero  thickness. 

37.  The  properties  of  interfftce  dements  are  defined  by  an  interfile  normd  stiffiiess 
kn  and  an  inteifoce  shear  stiffiiess  kg.  These  values  of  stiffiiess  relate  the  average 
rdative  di^lacement  normd  to  the  interfoce  dement  ^  and  average  rdative  shear 
di^lacement  to  the  corre^nding  notmd  stress  Off  and  shear  stress  r  the  equa¬ 
tions 


=  (9) 

and 

r=*A  (JO) 

The  units  of  and  kg  are  force  per  cubic  length. 
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38.  The  initial  value  of  u  set  equal  to  1  x  10*  within  the  inrogram.  This  value  for 
kff  ensures  that  the  normal  relative  di^lacement  of  the  interftce  element  is  insignificant 
viHhen  English  units  (fiset,  pounds)  or  SI  units  (meters,  Idlonewtons)  are  used.  If  other 
units  are  used,  the  ^ue  of  the  nomud  stififiwss  may  need  to  be  duuiged  to  a  higher 
value. 

39.  Two  types  of  interfiice  shear  response  are  modded,  a  Inlinear  shear  stress-dis> 
placemoit  rdationship  shown  in  Figure  4a,  and  a  hyperbolic  shear  stress-diq>lacement 
relationship  shown  in  Figure  4b.  In  the  bilinear  model,  the  value  assigned  to  is  a 
constant  so  long  as  the  average  shear  stress  t  almig  the  interfoce  is  less  than  the  shear 
strength.  If  the  shear  strength  of  the  interfile  dement  is  fiilly  mobilized,  udiich  occurs 
\dien  t  is  equal  to  y,  is  s^  equal  to  zero.  When  the  normal  stress  0|,  is  greato’  than 
or  equal  to  zero,  the  value  of  y  is  given  by  the  relationship 

Tf=c^+a,tUiS  (11) 


where 

C|  -  cohesion  intocept  dong  the  interfoce 
5  «  angle  of  intemd  fiiction  dong  the  interfiuMi 
and  shown  in  Figure  Sa.  When  is  less  thar.  zero,  y  is  conqnited  using 


i^iere  c/  *=  tensile  strength. 

40.  Direct  shear  test  results  on  soU-to-concrHe  interfiice  and  soil-to  sted  interfiice 
by  Potyondy  (1961),  Clough  and  Duncan  (1S>69),  and  Peterson  et  d.  (1976)  have 
shown  that  the  value  of  5  is  proportiond  to  the  angle  of  iittemd  fiiction  of  the  soil. 
The  value  of  the  constant  of  proportionality  is  dependent  upon  both  the  type  of  soil 
and  the  type  of  materid  comprising  the  sur&ce  of  the  structure. 

41.  The  direct  shear  tests  performed  by  Clou^  and  Duncan  (1969)  and  Petown  et 
d.  (1976)  have  shown  that  for  some  materids,  such  as  sand-to<onCTete  iiit»fiu:es,  the 
iiiterfiu;e  re^nse  during  shear  is  nonlinear  and  dqrendent  upon  the  nomud  stress.  A 
nonlinear,  stress-dqMndent  hyperbolic  cur^  is  used  to  represent  the  rdationship  be¬ 
tween  shear  stress  and  average  relative  shear  displacemoit  devdoping  during  primary 
loading  of  the  interfiice  (Figure  Sa)  and  a  linear  shear  stress-rdative  displaconent  re¬ 
sponse  during  uidoading  or  rdoading  of  the  intCTfiwe.  The  stress-dqrendent  value  of 
kgf  is  computed  using  the  rdationship 
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h.  tfyperbotic  representation  of  lAe  ycartatkmqf  shear  stress  with  relative 
shear  di^lacement  (kg^  ^  tangent  intatfcKe  diecir  stiffness^ 
kgf  -  initial  intaface  diecr  stiffltess) 


Figure  4.  Bitinun  tend  hyperbolic  models  far  interface  dtear  stress^ktave  shear 
tSsfdacement  behavior  after  CUmgftandDimcan(1969) 
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.St,)' 


(13) 


where 

^  initial  interfile  shear  stififiiess 
>  fiuhire  ratio 
SLj »  stress  levd 

The  initial  inter&ce  shear  stiffiiess  kgf  is  equal  to 


(14) 


where 

Kj  -  interfiice  nsoduhis  number 
Ym»  -  unit  weight  of  water 
n/ « interfiice  modulus  exponent 

42.  The  proportion  of  mobilized  shear  ^rength  for  each  interfoce  dement  is 
reflected  in  the  vdue  of  the  stress  levd  SLf.  SL^  is  equal  to  the  current  shear  stress  t 
divided  by  the  stress  at  foihire,  y 


(15) 

Tr 

Tj-  is  computed  using  either  Equation  1 1  12.  SLf  ranges  in  value  between  zero  and 

one. 

43.  The  fiulure  ratio  rdates  the  ultimate  shear  stress  to  the  shear  stress  at 
fiuhire 


The  ultimate  shear  stress  is  the  asymptote  to  the  shear  stress-relative  shear  duq>iace- 
ment  hyperbola,  as  shown  in  Kgure  4.b.  Direct  shear  stress  tests  on  sand-to-concrete 
intofoces  by  Peterson  et  al.  (1976)  have  shown  the  value  of  ^caDy  ranges  in 
value  from  0.3  to  1.0. 
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44.  The  relationship  between  the  average  normal  stress  along  the  intofiux  aiul  the 
tensile  strength  is  shown  in  Figure  5b.  The  value  of  kff  is  a  constant  value  equal  to  1  x 
10^  when  a„  is  greater  than  or  equal  to  o/.  If  a„  is  less  than  o/,  k„  is  set  equal  to  zero, 
assuring  that  additional  tensile  stresses  do  not  accrue  upon  subsequent  loadings.  This 
procedure  allows  for  separation  to  occur  between  two  adjacent  regions  of  along 
interface  elements  during  the  course  of  an  incremmtal  analysis.  If  the  separation  closes 
during  subsequent  loadings,  k„  is  reset  to  1  x  10^. 


Base  Separation  Analysis  >  ALPHA  Method 

45.  The  alpha  method  is  used  to  analyze  walls  that  are  loaded  so  heavily 
horizontally  that  gaps  develop  along  the  inter&ce  between  the  base  and  the  foundation 
(Ebeling,  Duncan,  and  Clough,  1990  and  Ebeling,  Clough,  Duncan  and  Brandon 
1992).  In  the  finite  element  model  of  the  soil-structure  interaction  problems,  the 
boumlary  between  the  base  of  the  monolith  and  its  foundation  is  modeled  usif^ 
interfiu:e  elements.  This  version  of  SOILSTRUCT  checks  for  separation  of  the  base  of 
the  wall  fi'om  its  foundation  after  each  increment  of  loading.  Base  sq)aration  occurs 
when  tensile  stresses  within  an  interfiice  element  center  are  computed  to  exceed  the 
tensile  strength  along  the  interfiice.  At  this  stage  in  the  analysis,  the  alpha  method  is 
initiated  by  the  computer  code.  The  purpose  of  die  alpha  me^od  is  to  rklistiibute  the 
tensile  stress  and  shear  stress  along  the  separated  interfiice  element  to  the  adjacent 
dements.  Details  reading  this  numerical  method  are  summarized  in  the  following 
paragraphs. 

46.  After  each  load  step  is  completed,  all  interface  elements  are  checked  for  the 

development  of  tensile  stress  at  thdr  centers.  This  occurs  prior  to  updating  the  stress, 
displacement,  and  modulus  arrays.  If  no  tensile  stresses  are  found,  the  analysis 
proceeds  with  the  array  updates  and  the  next  increment  of  load  as  u»ial.  When  toisile 
stresses  are  observed  in  the  interfiice  dement(s),  the  following  series  of  computations 
(idealized  in  Figure  6  for  a  single  interfiu:e  donent  developing  tension  after  application 
of  the  fiiU  load  increment  ^  updates  or  the 

analyze  for  the  next  load  increment: 

g.  For  each  interfiu:e  dement  that  devdops  tentile  stress  at  its  center,  the 
fi^ction  of  the  iqiplied  force  that  would  result  in  zero  normal  center  stress 
is  computed.  This  fiiiction  is  referred  to  as  alpha  and  is  computed  for  all 
elements  that  devdop  tensile  stress  at  thdr  centos.  Each  of  these  alpha 
values  may  be  described  as  the  fiaction  of  the  total  force  which  results  in 
zero  overshoot  normal  fisree  for  that  interfiice  donoit.  The  linear 
relationship  is  due  to  the  incremental  linear  dastic  anal^s  procedure  used 
in  the  program.  Due  to  the  mechanics  of  the  crack  development,  the 
interfiu^e  element  dosest  to  the  heel,  vdiich  has  not  devdoped  tensile  stress 
at  its  cento-  in  a  previous  load  incronent,  possesses  the  smallest  value  of 
alpha. 
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a.  Strength  criteria 


b.  Normal  stress-normal  relative  displacement  behavior 


Figure  5.  Interface  strer^crimiaandrrormalstress-nonnal 
relative  displacement  relationdtip. 
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Alpha  Mthod 

(a)  Apply  full  load  increment  {Qq}  to  devekp  tension  in  single  intetface  element. 


(b)  Apply  force  loads:  a^Qo  +  02^0 


Figure  6.  Idealization  of  the  a^dta  method. 
(Ebeling,  Clou^  Duncan  and  Brandon,  1992) 


b.  The  load  vector  {Q)  is  thoi  redefined  as  the  product  of  the  innallest  value 
of  alpha,  ax,  and  the  original  load  vector  for  this  load  stq)  (Qo)- 

{2}  =  «,{&}  (17) 

A  new  analysis  is  conducted  using  the  reduced  load  vector,  and  all 
displacements  and  stresses  are  updated. 

£.  Zero  noimal  and  shear  stifi&iess  is  then  assigned  to  the  interfiu:e  dement. 
This  prevents  the  accrual  of  stress  described  in  step  d. 

d-  For  the  critical  inter&ce  element,  the  tendon  stress  at  its  center  is  now  zero, 
but  artificial  restraining  forces  are  needed  to  reduce  the  built-in  slwar  stress 
to  zero.  The  artifidal  restraining  force  Qg  is  computed  for  the  donent. 
Using  the  unbalanced  load  procedure,  Qg  is: 

(18) 

where  {or^}  is  the  element  shear  stress  to  be  restrained  by  the  nodal  forces. 
It  is  equd  to  the  actual  linear  shear  stress  distribution  across  the  dement. 

is  the  transformation  matrix  that  relates  demem  strains  to  the  nodal 
point  displacements.  In  the  program  code,  this  is  accomplished  by 
subtracting  this  equivalem  shear  stress  from  the  stress  regime  existing 
across  the  entire  elernem. 

fi.  Since  the  force  Qg  does  not  actually  exist,  an  analysis  is  made  of  the  entire 
mesh  to  these  nodal  poim  forces  applied  in  the  opposite  direction  at  the 
nodes.  All  displacements  and  stresses  are  updated.  Essentially,  the  shear 
stresses  are  ralistributed.  This  step,  applied  in  conjunction  with  step  d, 
maintains  equilibrium  within  the  ^stem. 

f.  The  next  firaction  of  the  initial  load  to  be  applied  is  computed.  The  load 
vector  {Q}  is  defined  as; 


{e}=«,(i-z®){a}  (19) 

where 

aj  =  the  smallest  of  ..le  ronaining  values  of  alpha 
£a  =  the  sum  of  pre>dousfy  tq)plied  alpha  values 
{Qq}=  tlw  smallest  of  the  remaining  values  of  alpha 
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The  analysis  is  then  conducted  using  {Q},  and  all  resulting  diqilacefnents 
and  stresses  are  updated. 

g.  Steps  c  through  f  are  repeated  until  the  total  original  load  vector  for  this 
load  step  {Qq)  is  ^>plied.  Then  the  conventional  analysis  is  resumed  with 
the  iq)plication  of  the  next  load  increment. 

For  additional  details  regarding  the  alpha  method  and  the  unbalanced  load  procedure, 
see  Part  IV  of  Ebeling,  Clough,  Duncan  and  Brandon,  1992. 

47.  A  bilinear  shear  stress-displacement  rdationship  shown  in  Figure  4a  arul  a 
bilinear  normal  stress-displacement  are  used  for  the  imerface  elonents  along  the  base 
of  the  wall.  Use  of  the  alpha  method  is  restricted  to  walls  with  horizontal  wall  to 
foundation  interfaces.  The  data  input  for  the  base  separation  analysis  u^g  the  alpha 
method  is  described  in  Section  4  of  Appendix  A.  These  interftce  elements  must  be 
assigned  consecutive  numbers  along  the  base  of  the  structure. 
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PART  IV:  MODELING  THE  FAR-FIELD  BOUNDARY 


Introdnctioii 

48.  This  version  of  SOILSTRUCT  has  the  capability  of  taking  into  consideration 
the  effects  of  the  &r-fidd  domain  by  using  boundary  dements  based  on  a  half-plane 
fundamental  solution.  Figure  7  illustrates  a  typical  soil-structure  interaction  problem 
where  &r-fidd  domain  is  modded  using  boundary  dements.,  vdiile  the  structure  and 
adjacent  soil  domain  is  modded  using  finite  dements. 

49.  This  sections  describes  the  procedure  used  to  modd  the  far-field  of  a  soil- 
structure  interaction  problem  udng  boundary  donents.  A  brief  description  of  the 
boundary  element  us^  is  discussed  in  the  next  section.  The  substructure  method  for 
coupling  the  boundary  element  and  finite  dement  system  is  described  in  succeeding 
sections. 


Boundary  Elements 

SO.The  boundary  element  method  is  a  numoical  technique  for  solving  boundary 
value  problems.  For  an  elasticity  problem,  the  method  involves  discretizing  the 
boundary  of  the  problem  into  a  series  of  elonents.  Over  each  of  these  elements,  the 
displacements  and  tractions  are  chosen  to  be  piecewise  irv  :  rpolated  between  element 
nodal  points.  An  integral  equation  is  obtained  from  the  go  timing  differential  equations 
by  successive  application  of  integration  by  parts.  This  i^Tt^gral  equation  is  applied  to 
each  nodal  point  on  the  boundary  and  tlw  resulting  integrals  are  evaluated  (usually 
using  a  numerical  quadrature)  over  each  element.  A  system  of  linear  equations  in  terms 
of  the  nodal  boundary  displacements  and  tractions  results.  Prescribed  boundary 
conditions  are  then  iqrplied.  The  resulting  ^stem  of  equations  is  then  solved  for  the 
unknown  nodal  boundary  displacements  and  tractions. 

til.  SOILSTRUCT  uses  the  direct  boundary  element  method.  This  method  results 
in  a  system  of  equations  of  the  form: 

[«]{»}  =  [G]{p}  (20) 

>^ere  {r/}  are  the  nodal  boundary  displacements  and  {p}  are  the  nodal  boundary 
tractions.liie  matrices  [H\  and  [G]  are  cdled  the  matrices  of  influence  coefficients  and 
depend  on  the  type  of fundament  solution  used  as  well  as  ^e  manner  in  \^ch  the 
displacements  and  tractions  vary  between  nodal  points.  A  more  detailed  discussion  of 
the  direct  boundary  element  me^od  can  be  obtained  from  Brdrbia  (1984). 

52.  The  boundary  element  formulation  used  for  SOILSTRUCT  is  based  on  the 
infinite  half-plane  fundamental  solution.  This  fundamental  solution  gives  the 
displacements  and  stresses  that  result  due  to  the  application  of  a  concentrated  load  to 
the  interior  of  an  elastic  half-plane.  The  solution  to  this  problem  was  first  presented  by 


25 


Mdan  (1932).  Howevw,  in  the  original  work,  there  seem  to  be  errors  which  were 
pointed  out  and  corrected  by  Telles  and  Brd^ia  (1981). 

53.  In  this  version  of  SOILSTRUCT,  linear  isoparametric  boundary  demoits  are 
used.  The  geometry  of  the  dement,  as  shown  in  Figure  8,  is  defined  using  two  points 
(xij'i)  and  (X2,>'2)-  Within  the  elonent,  the  cartesian  coordinates  (x.y)  interpolated 
fi'om  the  locd  coordinate  ^  using  the  relationship 


X 


isl 


where  m/i(^)  are  the  lagrangian  interpolating  fimctions  given  by 


(21) 


V'l  = 


(1-^) 

2 


(22) 


In  a  similar  manner,  the  displacements  and  tractions  along  the  boundary  are 
interpolated  in  terms  of  the  nodal  values  uring  the  relationship 


(23) 

(24) 

is|  is| 


vdiere  and  Uyj  are  the  x  and  y  componoits  of  the  displacement  at  node  /.  Similariy, 
tjgj  and  tyj  are  the  x  andy  components  of  the  traction  at  node  r. 

54.  One  limitation  of  the  boundary  dement  method  is  that  for  problems  involving 
non-homogenous  or  non-linear  matoial  properties,  the  solution  procedure  becomes 
extremely  complicated.  To  keep  the  solution  procedure  simple,  the  &r  field  is  assumed 
to  be  a  homogenous  linearly  elastic  (isotropic)  half-plane.  A  single  value  of  E  and  v  is 
specified  for  the  entire  far  field. 


Coupled  Boundary  Element-Finite  Element  Method 

55.  To  couple  the  boundary  elernem  method  with  the  finite  element  method,  the 
entire  fin  field  is  treated  as  a  substructure.  This  is  done  by  converting  system  resulting 
fi'om  the  direct  boundary  element  method  as  described  in  equation  (20)  into  the  form 
of  the  direct  stififiiess  method; 
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FAR  FIELD  DOMAIN 


Figure  7.  A  typiad  smlstructure  interaction  problem  with  die 
fear  field-domain  modeled  usiiig  boundary  elements. 


Figure  8.  Linear  boundary  element  in  cartersUm  and  local 

coordinates. 


] 


[^K''}=w 


(25) 


{»}  are  the  nodal  displacements,  {F}  are  the  nodal  forces,  and  [i^]  is  the  equivalent 
stiffiiess  matrix  for  the  boundary  element  system.  In  this  form,  the  boundary  element 
stifBiess  matrix  can  be  assembled  into  the  ^obal  stiffiiess  matrix  just  like  ar^  finite 
element. 


S6.  The  equivalent  stiffiiess  matrix  for  the  boundary  dement  system  is  givoi  by 

[ie]  =  [A/lGr'[/(]  (26) 

where  [G]  and  [H\  are  the  matrices  of  influence  coefficients  described  in  equation 
(20),  and  is  a  matrix  used  to  convert  nodal  tractions  to  nodal  forces  using  the 
following  relationship 

{f}  =  [A/]0>}  (27) 


57.  For  a  single  boundary  element  with  linearly  varying  tractions,  the 
transformation  matrix  [A^  is  given  by 


2  0  10 
0  2  0  1 
10  2  0 
0  10  2 


(28) 


^ere  /  =  •y/(xj  -x,)^  +(>2  -y,)* .  TheAf  matrix  for  the  entire  ^stem  may  be  obtained 
by  assonbling  together  Af  matrices  of  the  individual  boundary  elements. 

58.  SOILSTRUCT  is  an  incremental,  equivalent  linear  analysis  program.  The  mesh 
used  must  have  a  sufficient  number  of  dements  so  that  two  dimensional  finite  dements 
are  used  to  modd  the  r^ion  in  which  noidinear  stress-strain  re^nse  is  to  be 
antidpated  while  boundary  dements  are  used  to  modd  the  r^on  idim  linear  stress- 
strain  re^nse  is  to  be  antidpated. 

59.  The  stiffiiess  matrix  is  always  un^mmetric.  To  preserve  the  ^mmetty  of 
the  finite  donent  ^em,  only  the  symmetric  part  of  [J^]  is  used  in  SOILSTRUCT. 
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The  (i,/)ch  dement  of  the  ^nunetiic  sdffiiett  matrix  [J^f  can  be  obtained  uan^  the 
follow^  fdatkmahtp: 

i, '  =  !(*,+*,)  (29) 

Brebbia  and  Walker  (1972)  have  shown  that  the  symmetric  matrix  [i^]^  resuhiqg  from 
equadmi  (29)  ctMreqxmds  to  a  least  squares  approximation  to  the  unsymmetric  matrix 

A 

[K]. 

60.  Zarco  (1993)  pmnted  out  that  fi>r  a  boundary  element  fomnilation  using  the 

A 

Mdan  fiindamental  solution,  the  dements  of  the  skew-symmetric  part  of  {K}  vdien 
normalized  with  re4>ect  to  the  dia^mal  terms  of  [JH  in  general  are  very  small: 

50.3%  Vi,>  (30) 

For  a  proUem  with  symmetric  geometry  and  loading,  if  only  half  of  the  problem's 

A 

domain  is  modded  using  the  "Method  of  Images”,  the  resulting  stiffness  matrix  [K] 
beccMnes  more  un^mmetric : 

pC^5l.0%  Vr,y  (31) 

I  I 

However,  the  additional  asymmetry  introduced  by  the  "Method  of  Images”  does  not 

have  a  very  significant  effect  when  [iH  is  assetrd>led  into  the  global  finite  demem 
stiffiiess  matrix  [ilQ.  This  is  due  to  the  symmetric  positive  definite  property  of  the  finite 

elemem  stififiMSS  matrix.  If  the  unsymmetric  stifiBiess  matrix  [JT]  were  assembled  in  to 
the  global  finite  demem  stiffness  matrix,  the  skew-^mmetric  terms  of  the  resulting 
stiffiiess  nutrix  [A]  would  be  very  snuD: 

50.01%  Vf,y  (32) 

^Miere  is  the  (iV)th  dernem  of  the  global  stiffiiess  matrix  [A]  resulting  fi^m  the 
assembly  of  [A]  into  the  ^obal  finite  dernem  stiffiiess  matrix. 
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61.  Zarco  (1993)  devdqied  another  couf^ed  BEM-FEM  program,  BEFEC,  wdiich 

uses  the  full  unsymmetric  stiffiiess  matrix  [JC].  This  program  is  c^Mble  of  solving  non< 
associated  plasticity  problems  using  the  Newton-Raphson  method.  For  non-assodated 

plasticity  problems,  it  is  necessary  to  retain  the  unsymmetric  portion  [K]  of  because 
the  global  stiffiiess  matrix  [il^  resulting  from  the  finite  dement  method  is  also 
unsymmetric.  For  the  linear  case,  aiulyses  performed  by  Zarco  using  BEFEC  indicate 

discardirig  the  skew-symmetric  part  of  [J^]  has  very  little  effect  on  the  resulting 
displacements  and  stresses. 

62.  SOILSTRUCT  modds  the  nonlinear  stress-strain  behavior  using  a  series  of 
equivalent  linear  analyses  in  which  the  finite  dement  stiffiiess  matrix  is  symmetric. 

Only  the  symmetric  part  of  is  used  in  SOILSTRUCT  so  as  to  preserve  the 
symmetry  of  the  finite  element  system.  Because  the  arudyses  are  linear,  discarding  the 

skew-symmetric  part  of  [J0]  does  not  significantly  affect  the  resulting  displacement 
stresses. 
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PART  V:  USE  OF  SOOSTRUCT  PROGRAM 


Sign  C«av«atiM  aad  CMrdfauite  SjrMtai 

63.  All  input  data  and  results  are  wptdStd  using  a  ri^-4iaad  coordinate  qfstem;  die 
x-axis  being  hcmzoittal  and  positive  to  the  and  the  y-axis  being  vettkal  and  posi¬ 
tive  iqiwards.  The  sign  convention  for  stresses  acting  at  the  center  of  a  two-dimen¬ 
sional  element  is  shown  in  figure  9a.  Compressive  stresses  are  taken  to  be  positive. 

64.  Stresses  for  inter&ce  dements  are  defined  with  reflect  to  their  local  axes  along 
the  length  of  the  inter&ce  x*  as  defined  by  the  1  and  J  nodes,  and  normal  to  the  ^ment 
y*  as  shown  in  Rgure  9b.  Positive  nonnal  stresses  are  congiressive.  Positive  shear 
stresses  act  in  the  positive  x*  direction  akx^  the  length  of  the  inteifoce  as  shown  in 
IHgure9b. 

65.  Positive  fivces  are  taken  to  be  conyressive  in  all  one-dimensimial  bar  elements 
with  the  exceptitMi  of  bar  reinfbrcement  dements.  For  diese  dements,  positive  ftnices 
are  taken  to  be  tensile. 


Units 

66.  Any  consistent  set  of  units  can  be  used  with  SOfiLSTRUCT,  with  om 
cautionary  note.  The  nonnal  stiflfhess  of  inteifoce  elements  is  arbitrarily  set  to  a  value 
of  1  X  10",  independent  of  units,  as  discussed  in  paragraph  38. 

Capacity 

67.  The  oyadty  of  the  program  is  determined  by  die  size  of  variables  in  the 
common  block  and  dimension  statements.  The  ^obd  sdffoess  matrK  [K]  is  stored 
within  the  inogram  vector  SN.  The  tfimendon  SN  is  equd  to  1,600,000,  and  the 
maxmnim  nuihber  of  either  nodes  or  two-dimensiond  cfcments  is  3,000.  the  dimendmi 
fin-  all  of  the  arrays  are  given  throughout  the  input  guide.  Double  precision  is  specified 
fm*  an  arrays. 


68.  The  giobd  sdfifiiess  matrix  is  assemUed  in  vector  fixm  usyg  a  modified  sl^ine 
procedure.  Therefiire,  there  is  no  direct  cmrdation  between  die  size  of  SN  and  the 
maximum  mmfoer  of  nodes  or  dements.  If  the  requred  size  of  SN  exceeds  1,600,000 
during  execution  of  SOILSTRUCT,  the  program  wiD  temmiate  execution  and  prim  the 
size  required  fisr  the  SN  vector.  Ifthis  occurs,  two  stateoMntsimist  be  modified  within 
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a  Coordinate  system  and  sign  convention  for  posithfe  stresses 
for  two-dimensiorial  elements 


b.  DeftniHon  of  mterface  hcai  ases  and  dgft 
cotwention  for  positive  shear  stress 


Figure  9  Sign  cortventions  arid  coordinate  axesfor  stress  and  displacenients 

(CkMgh  1984) 


the  computer  code;  (a)  the  dimensioa  of  the  SN  vector,  fixind  in  the  main  program, 
must  be  increased  to  ^  required  value,  and  (b)  a  caD  stater.ent,  comparing  tte  re¬ 
quired  size  of  the  SN  vecUM’  for  the  problem  bdng  analyzed  to  the  actual  dimenskm  of 
must  be  revised  in  subroutine  DETNA. 
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APFENDK  A:  USER’S  GUIDE  FOR  PROGRAM  SOILSTRUCT 


1.  IDENUFICATION  card  format  (20A4) 

CoImbmi  Variable _ Explanation _ 

I- 80  KED  Analysis  identification 

2.  DATA  CONTROL  CARD  FORMAT(14I5) 

(a)  MESH  PARAMETERS:  AU  nodes,  two-dimennonal,  interface  and  boundary 
elements  to  be  used  in  the  analyas  must  initially  be  included  in  the  mesh;  additions  or 
deletions  are  not  allowed,  but  the  material  parameters  can  be  changed  to  make  them 
inactive.  One-dimensional  bar  elements  usually  are  added  in  subsequent  construction  stqjs, 
but  may  initially  be  included  in  the  mesh. 

Column  Variable _ Explanation _ 

1-5  NUMNP  Number  of  nodal  points  (3000  maximum). 

6-10  NUMEL  Number  of  demoits  initially  in  the  mesh,  excluding  bar, 

boundary  and  beam  elements,  but  including  interfiu:e  elements. 
NUMEL,  then  includes  NUMJT  (3000  maximum).  Interfiu:e 
elements  should  be  numbered  first. 

II- 15  NUMJT  Number  of  inter&ce  elements  (200  maximum) 

16-20  NUMBAR  Number  of  bar  elements  initially  in  the  mesh  (60,  including 

elements  assessed  in  subsequent  construction  steps). 

21-25  NBEM  Number  of  boundary  elements  (200  maximum) 

(b)  ANALYSIS  PARAMETERS 

Column  Variable  • _ Explanation _ 

26-30  NC  Number  of  loading  construction  steps  (60  maximum) 

31-35  NMOD  Modulus  specification  code 

==  0  if  modulus  calculation  codes  input  with  loading 
information  card. 

==  1  if  modulus  calculation  codes  input  with  modulus 
calculation  card. 


Column  Variable _ E»phinotion _ 

36^0  INTT  Initial  stress  input  code 

=  0  if  external  input  from  cards  or  ti^,  included  in  input 
sequence.  Use  for  reading  initial  stresses  and  initial 
displacements  from  a  restart  file. 

=  1  if  internally  generated  from  gravity  tum-on  analysis. 

=  2  if  the  initial  stresses  and  displacements  are  to  be  set  equal 
to  zero  within  SOILSTRUCT 

41  >45  KI  Equal  to  zero  (int«&ce  is  automatically  activated).  KI  is  used 

when  INTT  is  not  equal  to  0.  If  INIT  =  0,  then  KI  can  be  left 

blank. 

46-50  IHORIZ  Ground  surfrice  inclination  code 

=  0  horizontal  ground  surface. 

Vertical  stresses  computed  from  gravity  tum-on.  Horizontal 
stresses  are  compute  assuming  Kq  =  v/(l-v).  Unless  Kq  is 
specified. 

=  1  sloping  ground  surface.  Vertical  and  horizontal  stresses  are 
calculated  from  gravity  tum-on  analysis  linear  elastic  response 
soil  -  i.e.,  Kq  =  v/(l-v). 

NOTE:  If  INTT  is  not  equal  to  0,  1  or  2,  initial  stresses  are  generated  assuming  a 
horizontal  ground  surface,  horizontal  water  table,  and  =  Kg  ol. 

(c)  OUTPUT  PARAMETERS 

Column  Variable  Explanation 

51-55  ITRD  Analysis  printout  code  - 

=  0  if  initial  stresses  are  to  be  printed. 

=  2  if  initial  stresses  are  not  to  be  printed. 

56-60  ILIST  Element  and  nodal  point  card  data  printout  code  - 

=  0  if  not  printed 
=  1  if  printed 

61-65  IPUNCH  Code  for  a  Restart  (tape  9)- 

==  0  Restart  file  not  created 
=  1  Restart  file  is  created 
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(d)  BASIC  PARAMETERS 


FORMAT(2D10^ 


Column  Varuibic _ Endanntion _ 

I- 10  GAMW  Unit  weight  of  water 

I I- 20  PATM  Atmo^hoic  pressure 

The  unit  weight  of  water  and  the  atmoq>heric  pressure  are  included  as  basic 
parameters.  Either  English  or  SI  units  can  be  used.  All  data  must  be  compatible  with  input 
coordinate,  pressure  and  material  property  parameters. 

3.  MATERIAL  ALLOCATION  CARD  FORMAT(F6IS) 

Column  Variable _ Eiplanation _ 

1-5  NUMMAT  Total  numbo*  of  material  types,  including  both  two- 

dimensioiud  soil  or  construction  material  types  and  inter&ce 
material  types.  (40  maximum) 

6-10  NUMSOL  Total  number  of  material  types  excluding  the  interface  material 

types.  Thus  (NUMAT-NUMSOL)  must  equal  the  number  of 
inter&ce  material  types. 

11-15  NATYP  Material  type  numbo’  asagned  to  2D  dements  having  the 

properties  of  air.  Usually,  elements  that  will  be  built  are 
initially  identified  as  air  elements. 

16-20  NATYP2  Material  type  number  assigned  to  interfiu^  dements  having  the 

properties  of  air.  Usually,  elements  that  will  be  built  are 
initially  identified  as  air  dements. 

21-25  NCTYP  Structural  material  type,  such  as  concrete  or  sheeqnling 

26-30  NBITYP  Backfill  material  type  1.  (refer  to  15b  on  fill  or  concrete 

placement) 

31-35  NB2TYP  BackfiU  material  type  2.  (refor  to  15b  on  fiU  or  concrete 

placement) 

4.  ALPHA  METHOD  PARAMTERS  FOR  BASE  SEPARATION  AND  FOR  2-D 
SOIL  ELEMENTS 

The  following  set  of  variables  control  the  modd  of  the  loss  of  contact  between  a 
structure  and  its  foundation  and  to  reduce  overshoot  stresses  in  2-D  soil  demoits  using 
the  alpha  method.  Note  that  only  horizontal  interfiice  dements  may  be  used  along  the  base 
of  the  structure  and  these  interfiice  dements  must  be  assigned  consecutive  numbers.  If  this 
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base  separation  routine  is  not  required,  assign  a  value  of  zero  to  NLOOP  and  proceed  to 
group  3  input  data. 

(a)  FIRST  CARD  FORMAT(2I540X,4I5) 

Column  Variable _ Eiplanation _ 

1-S  NLOOP  Maximum  number  of  iterations  using  alpha  method  for 

inter&ce  elemoits  for  each  load  step. 

6-10  KEYFRC  Linear  stress  increment  applied  to  inter&ce  elonents  during 

transfer  of  shear  stress  to  adjacent  elements 
=  0 

=  1  for  uniform 

31-35  NOELF  1  First  horizontal  inter&ce  element  checked  for  tension. 

36-40  NOELF2  Last  horizontal  intorfiu:e  element  checked  for  tention. 

41-45  NGROUP  Number  of  groups  of  interface  elements  for  which  resultant 

forces  are  to  be  computed. 

46-50  NLOOP2  Maximum  number  of  iterations  using  alpha  method  for  2-D 

soil  elements  for  each  loading  step. 

(if  the  alpha  method  for  2-D  soil  elements  is  not  required,  assign  a  value  of  zero  to 
NLOOP2) 

NOTE:  If  both  NLOOP  and  NLOOP2  are  non-zero,  then  the  smaller  of  the  two  will 
determine  the  maximum  number  of  iterations  using  the  alpha  method. 

Input  SECOND  and  THIRD  cards  for  each  plot  group,  NGROUP. 

(b)  SECOND  CARD  FORMAT(2I5) 

Column  Variable _ Explanation _ 

1  -5  NJELGRP  Number  of  elements  in  Plot  group. 

(40  maximum) 

6-10  NORIGIN  Node  about  vsdiich  Moments  will  be  summed. 

(c)  THIRD  CARD  FORMAT(16I5) 

Column  Variable _ Explanation _ 

1-5  NJEL  Element  Numbers  in  Plot  Group  NJELGRP.  (800  maximum) 

(NJELGRP)  Repeat  until  NJELGRP  elements  have  been  inputted. 
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4.1  CONTROL  PARAMETERS  USED  IN  CHECKING  EXCESSIVE 
INCREMENTAL  DISPLACEMENTS 

The  following  variables  control  the  routines  that  check  if  too-large  displacements  have 
occurred  in  the  mesh. 

(a)  FIRST  CARD  FORMAT(3I5^D10.14S) 

Column  VariaMe _ Eapiauation _ 

1-5  ICKCOL  Code  that  specifies  the  type  of  incremental  di^lacement  check 

to  be  performed. 

==  0  if  no  check  is  to  be  made 

=  1  if  a  check  is  to  be  made  for  excessive  x-displacements 
=  2  if  a  check  is  to  be  ntade  for  excessive  y-  displacements 
=  3  if  a  check  is  to  be  made  for  excessive  incremental  x  or  y- 
displacements 

6-10  LODCHK  Number  of  the  first  load  step  to  check  for  excessive 

incremental  displaconents.  The  same  check  will  be  performed 
in  all  subsequent  load  steps. 

11-15  NNDCHK  Total  number  of  nodal  points  to  be  checked  for  excessive 

incremental  displacements.  If  ICKCOLN),  the  value  of 
NNDCHK  will  not  be  used.  If  ICKCOL  does  not  equal  0  and 
NNDCHK=0,  dl  nodes  in  the  mesh  will  be  checked  for 
excessive  incremental  displacements. 

16-20  XTOL  Tolerance  value  for  nodal  incremental  x-displacement 

21-25  YTOL  Tolerance  value  for  nodal  incremental  y-displacement 

26-30  ICOLST  Stop  run  code 

=  0  Program  will  not  be  stopped 

=  1  Program  will  not  be  stopped  if  incremental  displacements 
of  specified  nodes  exceed  the  values  given  by  XTOL  and 
YTOL.  The  current  stresses  and  displacements  will  be  output 
prior  to  the  termination  of  the  analyris.  (The  results  will  be 
those  just  before  the  occurrence  of  excessive  displacements.) 

(b)  SECOND  CARD  FORMAT(16I5) 

Required  only  if  NNDCHK  does  not  equal  0.  The  nodal  point  numbers  of  nodes  (a 
total  of  NNDCHK)  to  be  checked  for  excesave  displacements  are  specified  in  the 


A5 


fi^wiQg  cards.  A  total  of  16  nodal  points  numbers  are  siq^died  on  one  card.  The  number 
can  be  input  in  any  order. 

Cduain  Variable _ Eiplanatioii _ _ 

1-S  NDCHK(I)  Code  that  specifies  the  type  of  incremental  di^laconent  check 

to  be  performed. 

=  0  if  no  check  is  to  be  made 

-  1  if  a  check  is  to  be  nude  for  excessive  x-displacements 
5.  LOADING  INFORMATION  CARD  FORMAT(6D^X45A4) 

One  card  is  supplied  for  each  loading  step.  One  to  three  loading/construction  modes 
can  be  included  in  each  loading  stq).  The  loading  or  construction  mode  codes  include; 

KCS(NCJ) _ DESCRIPTION _ 

1  Excavation  (equivalent  nodal  loads  can  be  applied  to  equal  incronents) 

2  Fill  placemoit  (subroutine  SOBSTP  cannot  be  used  in  conjunction  with 
the  ^  placemoit  procedures  of  subrmitine  BUILD) 

3  Seepage  loadir^  (equivaloit  nodal  loads  can  be  applied  in  equal 
increments) 

4  Deletion  of  bar  element  (force  in  the  element  can  be  iqiplied  in  equal 
increments) 

5  Installation  of  bar  element  (prestress  force  can  be  applied  in  equal 
increments) 

6  Boundary  pressure  loading  (equivalent  nodal  loads  can  be  applied  in  equal 
increments) 

7  Temperature  loading  (the  total  temperature  change  can  be  iqipUed  in 
equal  increments) 

8  Support  displacement  (the  total  di^lacement  can  be  applied  in  equal 
increments) 

9  Concoitrated  nodal  loads  (can  be  applied  in  equal  incremoits) 

10  Element  material  type  change 

As  indicated  in  the  listing,  iiq)ut  loads,  displacements  or  tenqjerature  changes  can  be 
analyzed  in  equal  increments,  or  substeps.  Subroutine  SUBSTP  generates  the  increments 
then  main  aiialyzes  all  increments  prior  to  analyzing  the  next  load  step.  With  one 


exception,  all  loading/constniction  modes  that  can  be  ^>plied  in  increments  or  substeps, 
can  also  be  qjplied  in  any  combination  in  any  load  step.  The  number  of  substeps,  however, 
will  be  the  same  for  all  loading  or  construction  m^es  included  in  the  load  stq).  The 
exception  is  temperature  loading;  if  a  temperature  change  is  ^)ecified,  and  a  given  number 
of  substeps  is  specified,  then  only  temperature  loading  can  be  specified  in  the  loading  stq) 
-  i.e.,  KCS(N,2)  and  KCS(N,3)  must  be  set  equal  to  zero,  or  left  blank.  If  the  number  of 
substeps,  NSBSP,  equal  to  zero,  then  temperature  loading  can  be  included  with  other 
loading/construction  modes  in  a  load  step. 


Since  the  same  input  format  is  used  and  similar  operation  are  performed  for  loading  or 
construction  modes  8  and  9,  the  following  rules  in  the  usage  of  these  two  codes  apply.  If 
orrly  concentrated  nodal  loads  are  specified,  use  mode  9.  If  only  support  displacements  are 
specified  use  mode  8.  If  both  loads  and  displacements  are  specified,  use  mode  8. 


Column _ Variable _ Expl«n*tion _ 

1  NALFPT(NC)  Printout  code  for  alpha  method. 

=  1  A  message  that  gives  details  on  the  alpha  method 
analyses  will  be  output  for  this  load  step.  Drafts  include 
which  element(s)  failed,  the  highest  stressed  element,  and 
the  value  of  alpha  for  this  element. 

2-3  KCS(NC,1)  First  loading/construction  mode  code 


T 


4-6 


KCS(NC,2)  Second  loading/construction  mode  code 


7-9 


KCS(NC,3)  Third  loading/construction  mode  code 


10-12  NUMIT(NC)  Number  of  iteration  for  the  loading  step.  NUMIT(NC) 

applies  to  each  substep  if  substeps  are  specified.  A 
maximum  of  10  iterations  can  be  spedfied.  NUMIT(NC)  = 
0  is  the  same  as  NUMIT(NC)  =  1 

13-15  NUMSS(NC)  Number  of  substeps.  A  maximum  of  10  substeps  can  be 

specified.  NIJMSS(NC)  can  not  be  zero. 


16-18  MOD(l,NC)  Modulus  calculation  code- 

=  1  if  a  loading  modulus  is  to  be  calculated 
=  2  if  an  unload  reload  modulus  is  to  be  calculated. 

=  0  if  the  computer  code  is  to  decide  the  type  of  modulus  to 
be  calculated.  In  tlus  case,  if  the  most  recently  calculated 
maximum  shear  stress  for  an  element  is  less  than  all 
previous  values  of  maximum  shear  stress,  an  unload-reload 
modulus  is  calculated.  Otherwise  a  primary  loading 
modulus  is  calculated 
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Input  is  required  here  only  if  NMOD  -  0.  AH  material  types  other  than  inter^  or 
elenienr*,  are  given  one  of  the  above  codes.  If  NMOD  ■  0  and  NC  •  0,  as  nnght  be 
case  fcMT  an  analysis  of  iiutial  stresses,  MOD(l,l)  is  set  equal  to  zero,  or 
SOILSTRUCT  code  decides. 

Column  Varinblc _ Eiplanatien _ _ 

19  IPRT(1,NC)  Print  code- 

~  0  print  displacements  and  stresses  for  fiiud  heraticHi  only 
■■  1  fffint  interfoce  stresses  for  each  iteration,  and  stresses 
for  all  dements  and  diq)laceiiients  for  final  iteration  only 
«  2  print  stresses  for  aU  dements  fin'  each  iteration,  and 
displacements  for  find  iteration  rmly 
«  3  ivint  stresses  and  diq)lacements  for  each  iteration 

20  1PLT(1,NC)  Code  for  creating  a  NISA  plot  files  (tapes  12  and  22) 

-ONISA  plot  file  not  created 
»  1  NISA  plot  file  created 
NISA  geometry  on  tape  12 
NISA  geometry  and  output  (Hi  tape  22 
21-80  HEDCS(NC)  Description  of  loading  stq> 

KCS(NC,1),  KCS(NC,2)  and  KCS(NC.3)  can  be  input  in  any  mimericd  order,  but  ^ 
modes  are  processed  in  ascending  mimericd  order:  If  second  and/or  third  loading 
construction  nixies  are  not  required  then  KCS(NC,2),  KCS(NC,3)  shcnild  be  set  etiud  to 
zero  or  left  blank. 

6.  MODULUS  SPECmCATION  CARD  FORMAT(4012) 

This  card  is  required  only  if  NMOD  *  1  and  NC  not  eijud  to  0.  A  card  is  repiired  for 
each  loading  step,  1  to  NC.  In  this  option,  values  of  the  mcxiijlus  qxxification  code  are 
qiecified  for  each  materid  type  (and  thus  each  dement,  exdudiiig  bar  and  interfile 
dements),  r^ardless  of  the  change  in  maxiimim  shear  strain  that  mmy  have  occurred. 

Column  Variable _ Eiplanation _ 

2,4,6,..  MOD(I,NC)  Modulus  cdculation  code  for  each  materid  type  (1  to 

NUMSOL)  for  the  first  load  step.  Separate  cards  are  required 
for  each  lotd  stqi.  Columns  not  used  can  be  left  blank. 

7.  MATERIAL  PROPERTY  CARDS 

These  cards  are  used  only  fi>r  two-dimensicMid  dements;  the  first  and  second  cards  must 
be  supplied  for  each  materid  type,  excludii%4Mur  and  interfile  materid  types.  The  first  and 
second  cards,  as  a  pair,  are  suppfied  in  order  of  materid  type  mindier  N  *  1  to  NUMSOL. 
Information  or  properties  not  required  for  a  materid  type  can  be  left  blank. 
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Hi 


(a)  FIRST  CARD 


FORMATXTDIO^IO) 


Colunia 

Variable 

Ezplanation 

1-10 

GUI(N) 

Poison's  ratio  before  fiuhire 

11-20 

GUF(N) 

Poison's  ratio  at  fiuhire  (no  greater  than  0.49) 

21-30 

GAM(N) 

Total  or  bouyant  unit  wdght  (always  qiecified  r^ardless  of 
drained  or  undrained  material  bdiavior) 

31-40 

41-50 

FR(N) 

AO(N) 

The  fiuhire  ratio 

Coeffident  of  lateral  earth  pressure  at  rest,  K,.  as  pertaining 
to  effective  stresses. 

51-60 

PHI(N) 

Friction  angle  in  d^rees 

61-70 

XXP(N) 

Expement  n  is  the  inhiaJ  tangent  and  unload-rdoad  modulus 
expressions.  Its  value  is  assumed  to  be  indqiendent  of  loadir^ 
m^e.  For  a  linear  elastic  material  n  must  be  0.  For  saturated 
soils  vriien  PHI(N)  «  0.,  N  should  normally  be  0.0001. 

71-80 

IDRAIN(N) 

Material  bdiavior  code  - 
»  0  if  undrained 
=  1  if  drained 

(b)  SECOND  CARD 

FORMAT(3D10J4>20.53D10.0) 

Column 

Variable 

Ezplanation 

1-10 

HCOEF(N) 

Coeffident  Km  in  the  initial  taitgent  modulus  eiqiression 

Ei  =  P^(a3/P,)n 

11-20 

ULCOEFF(N) 

Coeffident  Km  in  the  unload-rdoad  modulus  expression 

21-30 

COHE(N) 

Undrained  strength  or  cohesion 

31-50 

E(N) 

Taiigent  modulus  at  fiuhire  finr  isotropic  nonlinear  (non¬ 
elastic)  materials.  Youths  modulus  for  elastic  nuderials 

51-60 

ALPHACN) 

Coeffident  of  linear  thermal  eaqraiuion  for  structural 

dement;  zero  ot  blank  otherwise 
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Colwmn  Variable _ Exphumtion _ _ 

61-70  E1MIN(N)  Miiiiimiin  initial  tangent  modulus  for  isotropic  nonlimar 

(non-dasde)  mattfials.  Zero  or  left  blank  for  elastic  materials 

71-80  TENS(N)  Nfinimum  allowable  value  of  the  minor  prindple  stress  for 

isotropic  nonlinear  (non-elasdc)  nnaterials.  If  toisile,  input  as 
a  n^ative  value.  Zero  or  blank  for  dastic  materials. 

8.  INTERFACE  PROPERTY  CARDS 

One  pair  of  cards  are  supplied  for  each  intofue  material  type,  N  =  1  to  NUMJT.  If  no 
imerfiu:e  materials,  no  cards  are  required. 

(a)  FIRST  CARD  FORMAT(6D10.S4S) 


Column 

Variable 

Explanation 

1-10 

PHJ(N) 

Imerfiice  fiiction  an^e  in  d^rees. 

11-20 

RKS(N) 

Shear  stiffiiess  k,  before  &ilure 

21-30 

RKN(N) 

Normal  stifihess  before  ftilure 

31-40 

COJ(N) 

Inter&ce  cohesion 

41-50 

FRJ(N) 

Failure  ratio  Rfi,  FRJ  =  0.0  for  bilinear 

51-60 

TENSJ(N) 

Tendle  strength  of  interface  material  Of. 

61-65 

IADJMT(N) 

Material  number  of  2-D  soil  elonem  next  to  the  imerface 

(b)  SECOND  CARD 

FORMAT(4D10.5) 

Column 

Variable 

Explanation 

1-10 

STFSMN(N) 

Shear  stifihess  after  ftulure. 

11-20 

STFMN(N) 

Normal  stifl&iess  aft^  ftulure. 

21-30 

RKJ(N) 

Modulus  Number  for  interftice  dement  Kj. 

31-40 

XXPJ(N) 

Modulus  Exponent  for  imerfttce  dement  /r,. 

AlO 


9.  NODAL  POINT  CARDS 


F01IMAT(I10,4D10^) 


One  card  is  supplied  for  each  node.  The  min^bering  of  nodal  points  must  be  sequential 
and  some  of  the  nodes  can  be  omitted.  Those  nodes  omitted  are  automatically  goierated 
by  ttw  program  at  equal  spacings  between  those  q)edfied.  The  first  and  last  nodes  nuist 
always  be  specified.  Note  that  DP(N)  and  PP(N)  are  automatically  generated  in  equal 
increments  for  those  nodes  omitted  (see  Figure  Al) 


Column  Variable _ Eiplanation 


1-10 

N 

Nodal  point  ruimbo'. 

11-20 

X(N) 

X  coordinate,  positive  to  right 

21-30 

Y(N) 

Y  coordifude,  positive  upward 

31-40 

PP(N) 

Pore  pressure  in  head  of  water,  zero  or  blank  of  not  specified. 
Pore  pressure  must  not  be  specified  for  undrained  matoials, 
but  must  be  specified  for  drained  materials. 

41-50 

DP(N) 

Change  in  pore  pressure  in  hiead  of  water  for  soil  elements; 
change  in  temperature  for  linear  elastic  structural  material; 

zero  or  blank  otherwise. 

10.  BOUNDARY  CONDITION  CARDS 

Cards  1  through  8  are  supplied  as  required  to  specify  restraints  of  boundary  nodes.  If 
there  are  no  nodes  restrained  in  the  mode  specified  for  cards  2, 3,  or  4,  then  the  card  is  not 
required.  For  a  given  card,  specified  nodes  N=l,  to  NOY,  NOX,  or  NOXY  must  be  in 

sequential  order. 

(a)  FIRST  CARD 

FORMAT(3I5) 

Column 

Variable 

Explanation 

1-5 

NOY 

Number  of  nodal  points  fixed  against  Y-movement  only. 

6-10 

NOX 

Number  of  nodal  points  fixed  against  X-movement  only. 

11-15 

NOXY 

Number  of  nodal  points  fixed  against  both  X-  and  Y- 
movement. 

(b)  SECOND  CARD 

FORMAT(16I5) 

Column 

Variable 

Explanation 

I'S  IC(N)  Nodal  point  number  of  the  first  nodal  point  fixed  against  Y- 

movement 


All 


NODAL 


ELEMENT 

NUMBER 


PRIMED  NUMBERS 
INDICATE  NODE  OR 
ELEMENT  NUMBERS 
GENERATED  BY  THE 
COMPUTER,  USING 
THE  INPUT  (UNPRIMED) 
NODE  AND  ELEMENT 
NUMBERS. 


FigureAl.  Exanq^s  cf  number  nodes  and  eUnurtts 
if  generated  by  the  program  (Clough  1984) 
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Additional  nodal  points  fixed  against  Y-movemoit,  N  -  2  to  NOY,  are  q)edfied  in  the 
next  IS  five-cohunn  fidds  and  cm  addhicmal  cards  as  required. 

(c)  THIRD  CARD  FORMAT(16I5) 

Cdumn  Variable _ Explanation _ 

1-5  IC(N)  Nodal  point  number  of  the  first  nodal  point  fixed  against  X- 

movement 

Additional  nodal  points  fixed  against  X-movement,  N  »  2  to  NOX,  are  specified  in  the 
next  IS  five-column  fields  and  on  additional  cards  as  required. 

(d)  FOURTH  CARD  FORMAT(ldI5) 

Column  Variable _ Explanation _ 

I-S  IC(N)  Nodal  point  numb«‘  of  the  first  nodal  poim  fixed  against  both 

X-  and  Y-movement 

Additional  nodal  points  fixed  against  both  X-  and  Y-movonent,  N  =  2  to  NOXY,  are 
specified  in  the  next  IS  five-column  fields  and  on  additional  cards  as  required. 

11.  ELEMENT  CARD  FORMAT(6I5) 

Once  card  is  supplied  for  each  interface  or  two-dimensional  element;  bar  elements  are 
not  included  in  this  series  of  cards.  All  interfiu:e  elements  are  supplied  in  sequential  order 
first,  followed  by  two  dimensional  elements,  also  in  sequential  order.  (If  Build'  is  used; 
place  imerfiu:e  elements  which  will  be  buih  up  at  the  end  of  joint  Ust.  i.e.  start  numbering 
these  elements  with  N  =  NUMJT  -  NJTFIL).  Thus,  inter&ce  elements  must  be  numbered 
fi'om  N  =  1  to  NUMJT,  and  two-dimensional  dements  fi’om  N  =  (NUMJT+1)  to 
NUMEL. 

Nodal  point  numbers  must  be  specified  consecutivdy,  processing  coumer-clocicwise 
around  the  element  as  shown  in  Figure  A2  The  first  and  last  nodal  point  numbers  specified 
for  interfi^  elements  must  have  the  same  coordiiunes.  Triangular  two-dimensional 
elements  having  four  different  nodal  point  numbers  nuiy  not  be  used;  the  first  and  last 
point  numbers  of  a  triangular  elemoit  must  be  identical. 

Element  numbers  in  a  row  may  be  omitted,  in  whidi  case  the  omitted  dements  will  be 
goierated  by  incrementing  the  element  number  N  and  the  nodal  point  numbers  I,  J,  K,  and 
L  by  one,  and  by  assigning  the  same  material  type  number  as  specified  for  the  last  denmnt. 
The  first  and  last  dements  in  the  row  must  be  q>edfied.  If  no  dements  are  omitted,  the 
element  numbering  may  be  done  in  ai^  ordo*,  provided  all  into'fiue  dements  are 
numbered  first. 
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Etement  Numbo' 


VariaUe 


N 


n/N,  1)  Nuiri)er  of  Nodal  Poim  I. 

IL(N,2)  Number  of  Nodal  Point  J. 

IL(N,3)  Number  of  Nodal  Point  K. 

IL(N,4)  Nund)er  of  Nodal  Pcmt  L. 

nXN,S)  Material  Type  Number 


(a)  ACCEPTABLE 


THICKNESS' 


(b)  UNACCEPTABLE 


Figure  A2  Nodal  point  numbering  of  isopawneOic 
QMS  elements  and  interface  elements  (Clougfi  1984) 


12.  BAR  ELEMENT  CARD 


FORMAT(41S^D10.74D10.14>10.5) 


One  card  is  nipplied  for  each  bar  element  initially  in  the  mesh  or,  if  a  continuation 
analysis,  added  in  an  incremental  loading  step  of  the  previous  analysis.  Note  that  for  a 
continuation  analysis  this  card  is  not  automatically  generated.  Elements  are  numbo’ed 
sequentially  from  N  1  to  NMBAR 


Column  Variable _ Eiplanation 


1-5 

N 

Element  Number 

6-10 

IB(N,1) 

Number  of  nodal  point  I 

11-15 

IB(N,2) 

Number  of  nodal  point  J 

16-20 

IB(N,3) 

Spring  response  type  code- 
=  1  if  both  compression  and  tension  of  bar  allowed 
=  2  if  only  compression  allowed. 

~  3  if  only  tension  allowed. 

21-30 

BAR(N,!) 

cos  a 

31-40 

BAR(N,2) 

sin  a 

41-50 

BAR(N,3) 

Prestress  force  in  the  bar  element.  The  force  must  be  applied  as 
loads  at  nodes  I  and  J  using  loading/construction  mode  9  in  a 
loading  step. 

51-60 

BAR(N,4) 

Stiffiiess  of  bar  element.  This  usually  computed  as  AE/L,  but 
the  mesh  length  (distance  from  node  I  to  n^e  J)  need  not,  and 
usually  does  not,  correspond  to  the  actual  length. 

61-70 

BAR(N,5) 

Displacement  of  bar  element  at  activation.  This  allows  for  a 
specified  degree  of  slack  at  the  strut  connection;  the  bar  wiU 
deform  BAR(N,5)  prior  to  its  stiffiiess  being  activated. 

Bar  elements  can  function  as  either  anchors  or  strut  (spring)  supports.  The  required 
parameters  are  dependent  on  the  type  of  bar  element  specific. 

If  a  strut  support  is  specified,  nodal  point  J  is  a  node  fixed  against  x-  or  y-  movement, 
depending  on  the  orientation  of  the  strut  being  modeled.  For  program  storage  efficiency 
the  number  of  node  J  should  be  as  close  as  posable  to  the  numl^  of  node  I.  Nodal  point  I 
represents  the  point  of  connection  between  the  wall  and  the  actual  strut.  Nodal  points  I 
and  J,  then,  are  not  necessarily  physically  connected,  since  the  element  stiffiiess  is  input 
independently.  Nodal  point  J  allows  the  force  at  the  J  node  to  be  carried  into  the  ^stem  as 
a  reaction  at  a  fixed  node.  This  is  consistent  with  the  typical  mesh  representation  of  one 
half  of  a  symmetric  ^cavation.  The  value  of  sin  a  and  cos  a  are  specified  according  to  the 
sign  convention  shown  in  Figure  A3.  The  values  are  input  to  represent  the  line  of  action  of 
the  strut  support,  and  do  not  need  to  correspond  to  the  relative  positions  of  the  I  and  J 
nodes. 
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sin  a  Si- 
cos  a  =  - 


sin  as-i- 
cos  a  =  + 


Figure  A3.  Sign  conventim  and  definition  of  local  axes 
for  bar  element  (Clough  1984) 

If  an  anchor  is  specified,  nodal  points  1  and  J  physically  represent  two  ends  of  the 
anchor,  and  must  be  restrained  appropriately.  The  values  of  sin  a  and  cos  a  must 
correspond  to  the  relative  portion  of  the  I  and  J  nodes  representing  the  ends  of  the 
anchor.  Stiffiiess  is  computed  as  AE/L,  with  L  bong  the  distance  between  nodes  I  and  J 
and  either  A  or  £  alter^  to  give  the  correct  sfiffiiess.  StifBiess  for  an  anchor  or  a  stmt 
support  is  input  as  force  per  length  of  wall. 

For  either  element  type,  ^)edfying  the  prestress  force  does  not  rq)ply  the  force.  The 
concentrated  force  loading/constmction  m^e  must  be  used  for  this  purpose.  Thus,  bar 
elements  initiaUy  in  the  mesh  cannot  carry  prestress  force,  since  it  is  not  q}plied  by  a 
gravity  turn  on  aruilysis. 

13.  BOUNDARY  ELEMENT  CARD 

The  maximum  number  of  boundary  elements  that  can  be  specified  is  200.  Figure  A4 
shows  that  the  nodes  within  each  boundary  dement  are  numbered  so  that  normal  to  each 


A16 


towards  the  element  domain.  The  local  numboing  of  the  nodes  in  boundary  demon 
system  must  conform  to  this  confluency  convention  illustrated  in  Figure  AS. 

(a)  FIRST  CARD  (NODAL  COORDINATE  DATA)  FORMAT(16I5) 

One  card  is  supplied  to  specify  the  coordinates  of  the  global  nodal  point  numbers  in  the 
boundary  element  system.  The  specified  nodes  N==1,[NBEM+1]  must  be  in  sequential 
order. 

Column  Variable _ Explanation _ _ _ 

1-5  IX(N)  Global  nodal  point  number  of  first  node  in  boundary  element 

system. 

Additional  global  nodal  point  numbers,  N-2  to  NBEM,  are  specified  in  the  next  IS  five 
column  fields  and  on  additional  cards  as  required. 

(b)  SECOND  CARD  (SYMMETRY  CONDITIONS)  FORMAT(2I5) 

One  card  is  supplied  to  specify  the  ^mmetry  conditions  of  the  boundary  element 
system.Because  of  the  formulation  which  assumes  a  half-plane  with  the  surfiu^e  parallel  to 
the  x-axis  as  the  domain  for  the  boundary  element  system,  only  symmetry  about  the  y-axis 
is  considered. 

Column  Variable _ Explanation  _ 

1-S  NSYM  Number  of  nodes  in  boundary  element  system  lying  on  axis  of 

symmetry  (Y-axis).  If  NSYM  =  0,  no  symmetry  in  problem. 

(c)  THIRD  CARD  (SYMMETRY  COORDINATE  DATA)  FORMAT(16I5) 

One  card  is  supplied  to  specify  the  local  nodal  point  numbers  in  the  boundary  element 
system  which  lie  on  the  axis  of  symmetry.  If  NSYM  =  0,  skip  this  card.  The  specified 
nodes  N=1,NSYM  can  be  input  in  any  order. 

Column  Variable _ Explanation _ _ 

1-S  ISYM(N)  Local  nodal  poim  number  of  node  in  boundary  element  system 

lying  on  axis  of  symmetry  (Y-axis),  1  ^  ISYM(N)  ^ 
[NBEM+1] 

Additional  local  nodal  point  numbers  in  boundary  element  system  lying  on  the  axis  of 
symmetry  are  specified  in  the  next  IS  five-column  fields  and  on  additional  cards  as 
required. 
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INCORRECT  NUMBERING 
SCHBME 


Figure  A4.  Nodal  point  numbering  of  linear 
bouruiary  element. 


Figure  A5.  Local  rmmbering  of  bouruiary  element  nodal  points 
in  a  coupled  boundary-finite  element  ^stem. 
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(d)FOURTH  CARD  (MATERIAL  PROPERTIES)  FORMAT(2D103) 

One  card  is  supplied  to  q)ecify  the  material  properties  of  the  entire  boundary  dement 
^em. 

Column  Variable _ Explanation _ 

I- 10  E  Young's  Modulus  for  fiu*  fidd. 

II- 20  PR  Poison's  ratio  for  fiur  fidd. 

Since  only  one  set  of  E  and  v  must  be  specified  for  the  entire  fiur  fidd,  it  is  suggested 
that  the  average  value  of  E  and  v  at  the  bottom  of  the  mesh  be  used. 

14.  CONTINUATION  OR  INITIALIZATION  CARDS 

This  iiq)ut  is  i4)pUed  only  if  INTT  »  0;  it  can  be  supplied  fi’om  external  disk  storage  or 
punched  cards  input  fi'om  a  preceding  analysis.  Ii^ut  format  is  the  same  for  both  input 
modes.  This  option  is  provided  so  that  a  required  sequence  of  loading  stq>s  can  be 
stopped  at  an  intermediate  step,  thoi  restart^  fi-om  that  stq>  without  repeatii^  the 
complete  aiudysis.  These  cards  may  also  be  used  to  specify  values  of  particular  variables  in 
an  initial  atud^s  without  using  the  gravity  tum-on  procedure  followed  with  INTT  ==  1  or 
the  special  procedure.  Similarly,  vdues  assigned  to  spedfic  variables  can  be  changed  if  the 
sequence  of  loading  is  stopped,  prior  to  restart  analyds. 


(a)  FIRST  CARD  (ELEMENT  INFORMATION)  FORMAT(5IS) 


Column 

Variable 

Explanation 

1-5 

NUMEL 

Number  of  donents  in  the  mesh,  excluding  bar  elements,  but 
including  interfiice  dements. 

6-10 

NUMJT 

Number  of  inter&ce  donents. 

11-15 

NUMBAR 

Number  of  bar  donmits,  including  those  initially  in  the  mesh 
and  those  added  in  previous  loading  steps  (if  a  restart 
analysis). 

16-20 

NUMNP 

Number  of  nodal  points. 

21-25 

NSTART 

Interfiice  donent  code- 

=  0  if  all  intofiice  elements  are  active  (i.e.  no  shear  or  tension 
fiulures).  IFLOLD  and  IFAIL  keys  are  set  equal  to  zero  for  all 
interfile  elements  within  the  program.  See  card  groups  g  and 
h. 

==  1  if  IFLOLD  and  IFAIL  keys  are  to  be  read  for  all  interfiice 
elemoits  was  generated. 
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Whoi  NSTART  *  1,  the  keys  IFLOLD  and  IFAIL  have  in  a  previous  analysis  and  are 
included  as  standard  output  to  tape  9  when  the  disc  storage  option  (see  LOADING 
INFORMATION  CARD)  was  activated. 

(b)  SECOND  CARD  (STRESS  INFORMATION) 

INTERFACE  ELEMENTS  FORMAT(2(lP,lD13.64X)) 

2-DIM  ELEMENTS  FORMAT(4(lP4D13.64X)) 

Three  cards  are  supplied  for  each  interface  element,  in  numerical  sequence  N  =  1  to 
NUMJT.  For  each  inteiftce  element,  the  pair  of  normal  and  shear  stresses  at  node  I 
(SIGI),  the  center  (SIG),  and  node  J  (SIGJ)  are  supplied. 


Column 

Variable 

Explanation 

1-13 

SIG(N,1) 

Stress  in  the  x-direcdon  for  a  two-dimensional  element,  Ox, 
normal  stress  for  an  interface  element 

16-28 

SIG(N,2) 

Stress  in  the  y-direction  for  a  two-dimensional  element,  Oy, 
shear  stress  for  an  interfiu^e  element 

31-43 

SIG(N,3) 

x-y  shear  stress  for  a  two-dimensional  element,txy,  zero  or 
blank  for  interface  dement. 

46-58 

SIG(N,4) 

Maximum  previous  value  of  x-y  shear  stress  for  a  two- 
dimensional  elemoit,  zero  or  blank  for  interfiu:e  dement 

One  card  is  supplied  for  each  two-dimensional  element,  in  numerical  sequence  N  = 
(NUMJT +1)  to  NUMEL. 

(c)  THIRD  CARD  (NODAL  POINT  INFORMATION) 

FORMAT(3(lPlD14.74X) 

Information  for  three  nodal  points  is  supplied  on  each  card.  Nodal  points  are  specified 


in  numerical  order,  N  = 

1  toNUMNP 

Column 

Variable 

Explanation 

1-14 

DISPX(N) 

Displacement  in  x-direction 

16-29 

DISPY(N) 

Displacement  in  y-direction 

31-44 

PP(N) 

Pore  pressure  in  head  of  water 

(d)  FOURTH  CARD  (MATERIAL  TYPE  DESIGNATION)  FORMAT(15IS) 

Material  type  number  for  IS  interfiice  or  two-dimensional  elements  are  ^)ecified  on 
each  card,  hfoterial  type  numbers  for  elemoits  in  numerical  sequence,  N  =  1  to  NUMEL, 
are  specified. 
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Column  Varuible _ Eiplanation _ 

1-S  IL(^^,S)  Material  type  munber 

Material  type  number  for  the  next  14  demoits  are  supplied  in  the  imxt  14  five-cohimn 
fidds. 

Note  that  material  type  numbers  supplied  on  these  cards  supersedes  the  material  type 
numbers  specified  on  the  element  card  (section  1 1).  Thus  matoial  type  changes  can  be 
made  as  part  of  a  restart  analysis  rather  tluui  ii^uding  such  changes  in  a  loading  step  of  an 
analysis. 

(e)  FIFTH  CARD  (BAR  ELEMENT  INFORMATION) 

FORMAT(3I5»2D10.7^D10.14>10.5) 

This  card  is  supplied  only  if  bar  elements  are  included  (NUMBAR  >  0).  Information  for 
bar  elements  is  specified  on  each  card.  Information  is  supplied  for  bar  elements  in 
numerical  sequence,  N  -  1  to  NUMBAR. 


Column 

Variable 

Explanation 

1-5 

IB(N,1) 

Numbo’  of  nodal  point  I. 

6-10 

IB(N,2) 

Number  of  nodal  point  J. 

11-15 

IB(N,3) 

Spring  response  type  code  - 
=  1  if  both  compresdon  and  tension  of  bar  allowed 
=  2  if  only  compresdon  allowed 
=  3  if  only  tension  aUowed 

16-25 

BAR(N,1) 

cos  a 

26-35 

BAR(N,2) 

sina 

36-45 

BAR(N,3) 

Prestress  force  in  bar  element. 

46-55 

BAR(N,4) 

Stifihess  of  bar  donent,  (AE/L). 

56-65 

BAR(N,5) 

Di^lacemoit  of  bar  demem  at  activation. 

Note  that  these  parameters,  if  duuiged  for  a  restart  analysis  supersede  those  specified 
on  the  bar  demoit  card  (section  U).  Also,  if  bar  dements  are  inifially  included  in  the 
mesh,  and  an  inhialization  procedure  is  used,  then  tins  card  must  be  induded,  duplicating 
the  information  specified  in  the  Bar  Bonem  Card. 


FORMAT(SD10.4) 


(0  SECra  CARD  (INTERFACE  INFORMATION) 

This  card  is  supplied  only  if  interftce  denwnts  are  included  (NUMJT  >  0).  Infomiatum 
for  four  interfoce  dements  is  q)ecified  on  each  card.  InfiMmation  is  suf^lied  for  uiterfoce 
dements  in  numericd  sequence,  N  »  1  to  NUMJT.  If  no  interfiice  dements  are  present, 
proceed  to  group  12  input  data. 

Colmnn  Variable _ Eiplanation _ 

I- 10  STFS(N)  Shear  Stiffiiess  of  first  interfoce  dement 

I I- 20  STFN(N)  Normal  stiffiiess  of  first  interfiu:e  demoit 


Note  that  the  value  of  the  shear  stiffiiess,  if  changed  for  a  restart  analysis,  supersedes 
the  value  specified  on  the  Interfiu:e  Property  Card  (section  8).  Thus  the  interfiice  stiffiiess 
can  be  changed  as  part  of  a  restart  analysis. 

(g)  SEVENTH  CARD  (INTERFACE  KEY)  FORMAT(15I5) 

This  card  is  supplied  if  interfiice  dements  are  induded  (NUMJT  >  0)  and  NSTART  » 
1.  Information  for  fifteen  intofiice  donents  is  spedfied  on  each  card.  Information  is 
supplied  for  interfiu:e  elements  in  numerical  sequence,  N  -  1  to  NUMJT.  If  no  imerfiice 
elements  are  present  or  NSTART  =  0,  proceed  to  group  12  iiqiut  data 

Column  Variable _ Eiplanation _ _ 

1-5  IFLOLD(^0  Modulus  asdgnment  key - 

=  0  if  interfiice  dement  N  is  active 
-  1  if  tension  Mure  occurred 
=  2  if  shear  fiuhire  occurred 

(h)  EIGHTH  CARD  (INTERFACE  KEY)  FORMAT(15I5) 

This  card  is  supplied  only  if  NSTART  -  1  and  for  only  those  interfiu:e  elements  along 
the  horizontal  base  of  the  structure  used  to  modd  the  base  separation,  with  donoit 
numbers  from  NOELFl  to  NOELF2.  Information  for  fifteen  interfiice  elements  is  specified 
on  each  card. 

Column  Variable _ Eiplanation _ 

1-S  IFAIL(N)  Stress  transfer  k^- 

-  0  if  interfiice  dement  number  N  is  active 

-  1  if  transfisr  of  normal  and  shear  stress  in  fiuled  interfiice 
deriMnt  number  N  has  been  completed. 
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15.  LOADING  STEP  CARDS 


These  cards  are  required  only  if  NC  is  not  equal  to  0.  Cards  are  assembled  in  the  order 
specified  on  the  Loading  Information  Ciufi  (section  S):  cards  for  the  first 
loading/construction  mode  qiecified  for  the  first  loadiitg  stq),  and  so  on,  to  the  cards  for 
the  last  loading/construction  mode  q)ecified  for  the  last  construction  stq[>. 

For  a  given  loading  step,  the  lowest  numboed  loading/construction  mode  is  processed 
first,  but  the  analysis  of  the  loading  step  is  made  for  the  combined  effect  of  all 
loading/construction  imxles  included  in  the  loading  step.  Care  must  therefore  be  exo'cised 
in  specifying  some  loading/construction  modes,  such  as  material  type  changes  or 
concentrated  forces  representing  prestress  forces,  in  the  same  loading  stq)  with  other 
loading/construction  inodes. 

(a)  EXCAVATION 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  ==  1.  Input  is 
handled  by  subroutine  ^C AV.  Free  excavated  elements  and  common  excavated  elements 
are  input  separately.  A  fiee  excavated  dement  is  an  dement  specified  to  be  excavated  that 
has  sfi  boundaiy  in  common  with  demoit  specified  to  be  excavated  in  the  loading  step.  A 
common  excavated  dement,  therefore,  has  at  least  one  boundary  in  common  with  an 
unexcavated  dement  (Figure  A6) 

Interfiux  elements  can  only  be  included  as  firee  excavated  dements,  even  if  they  have  a 
boundary  in  common  with  an  unexcavated  dement.  Imerfiuxs  element  cannot  be  used  as 
interpolation  elements.  Free  excavated  doiwnts  (other  than  interpolation  dements)  can  be 
used  as  imerpolation  elements,  though  common  excavated  elements  are  more  commonly 
used. 

If  possible,  adjacent  common  excavated  dements  should  be  input  sequentially  as  this 
procedure  avoids  repetitive  computation;  nodal  loads  for  a  nodal  point  common  to  the 
two  sequential  elements  only  be  calculated  once. 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT(2I5) 

Column  Variable _ Explanation _ 

1-S  NFXEL  Number  of  fiee  excavated  dements. 

6-10  NXELCB  Number  of  common  excavated  demons. 
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Figuv  A6  Exanipk  excavatkm  had  Step  tkfina^frte- end 
cammon-excavat^  elements  in  relation  to  the  excavtOion 
bomday  (Ckm^  1984) 

(2)  SECOND  CARD  (FREE  EXCAVATED  ELEMENT  DATA) 

FORMAT(161S) 

Element  numbers  of  16  free  excavsted  elements  can  be  supplied  on  one  card.  A 
maximum  of  SO  can  be  q)ecified  in  on  knufing  step.  Element  numbers  of  all  free  excavated 
dements,  N  -  1  to  NFXEL  are  to  be  q)ecified. 

Colnain  Variable _ Explanation _ 

1-5  LNXEL(N)  Element  number  offirst  free  excavated  dement. 

Infrnniation  fr>r  the  next  IS  free  excavated  dements  is  supplied  in  the  next  IS  five- 
column  fidds. 

(3)  THIRD  CARD  (COMMON  EXCAVATED  ELEMENT  DATA) 

FORMAT(SIS) 

One  card  is  supfdied  fr)r  eadi  aunmon  excavated  dement,  N  *  1  TO  NXELCB.  A 
maximum  of  SO  conunon  excavated  dements  can  be  q>ecified  in  cm  loading  stq).  Loading 
codes  indude: 
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0  -  The  node  is  not  k>aded  by  e^vatkm  forces  and  is  not  common  to  both  an 
excavated  and  an  unexcavated  etement. 

1  -  The  node  is  loaded  by  excavation  forces  and  is  common  to  both  an  excavated 
and  an  unexcavated  dement. 

Note  that  I,  J,  K,  and  L  refo-  to  the  same  nodes  I,  J,  K,  and  L  specified  on  the  dement 
card  (section  1 1) 

Interpolation  dement  should  be  specified  in  a  criss-CToss  foshion  as  shown  in  Figure 
A7.  Further,  x-  and  y>coordinates  of  diagonal  dements  must  not  be  the  same.  If  these 
rules  are  not  adhered  to,  the  interpolation  routine  will  detect  singularity  and  processing 
will  stop. 


Column 

Variable 

Explanation 

1-5 

LUL(N,1) 

Element  number  of  the  first  common  excavated  elemoit.  This 
donoit  is  also  the  first  interpolation  dement. 

6-10 

LUL(N,2) 

Element  number  of  second  interpolation  dement. 

11-15 

LUL(N,3) 

Element  numbo*  of  third  interpolation  element. 

16-20 

LUL(N,4) 

Element  number  of  fourth  interpolation  element 

21-25 

LUL(N,5) 

Loading  code  for  node  I. 

26-30 

LUL(N,6) 

Loading  code  for  node  J. 

31-35 

LUL(N,7) 

Loading  code  for  node  K. 

36-40 

LUL(N,8) 

Loading  code  for  node  L. 

Loading  codes  are  specified  for  the  nodes  of  the  element  to  be  excavated  q)ecified  in 
columns  1>S. 

(b)  FILL  OR  CONCRETE  PLACEMENT 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  OR  KCS(N,3)  «  2.  Input  is 
handles  by  subroutine  BUILD.  The  types  of  dements  that  can  be  placed  include  structural, 
soil  and  inteifirce  elements.  Subroutine  SUBSTP  cannot  be  used  in  conjunction  with 
subroutine  BUILD. 
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niWT  MTERKIATION  EUaflENT 


m.  ABOVE  CONFIGURATIONS  b.  ABOVE  CONFIGURATIONS 

ARE  ACCEPTABLE  ARE  NOT  ACCEPTABLE 

Figure  A7.  lUustration of  mmbering of  interpoiation 
ekments for  calculation  of  stresses  at  excavation 
boundary  (Clougfi  1984) 

At  placement,  the  fill  denient(s)  is  assigned  a  low  modulus  and  the  surfiice 
diq>lacements  are  set  equal  to  zero.  Stresses  assigned  to  the  nei^  place  fill  are  based  on 
q,  =  K(Py,  where  Oy  is  equal  to  the  |Koduct  of  the  unit  wei^  of  the  fill  dement  times  the 
dqrth  bdow  the  sumce  to  the  center  of  the  dement. 
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(1)  FIRST  CARD  (CONTROL  DATA) 


FORMAT(5154>l<K2) 


Column 

Variable 

Explanation 

1-5 

NLEL 

Total  number  of  dements  to  be  placed  including  interftux  and 
structural  dements. 

6-10 

NJ 

Total  number  of  denents  to  be  placed  less  the  number  of 
inactive  interfiice  dements  to  be  placed. 

11-15 

NONP 

Number  of  nodal  points  within  the  placed  layer(s)  to  be 
assigned  zero  displacements.  This  include  all  nodal  points  of 
the  dements  to  be  placed  except  those  nodal  points  in  common 
with  an  already  existing  dement. 

16-20 

NCE 

Number  of  structural  dements  to  be  placed. 

21-25 

NBIE 

Number  of  backfill  type  1  dements  to  be  placed.  (The  numb«’ 
of  NB2TYP  elements  is  determined  internally.) 

26-35 

HTB 

New  y-coordinate  of  the  top  of  the  place  backfill. 

(2)  SECOND  CARD  (ELEMENT  NUMBERS)  FORMAT  (1615) 

Element  number  of  16  placed  dements  can  be  supplied  on  one  card.  A  maximum  of 
100  can  be  specified  in  one  loading  step.  Element  number  of  all  placed  elements,  N  =1  to 
NLEL,  are  to  be  supplied.  The  elements  must  be  read  in  by  material  type  according  to  the 
following  sequence: 

1.  structural  elements 

2.  backfill  type  1  elements 

3.  backfill  type  2  elonents 

4.  interfiice  elements  to  be  activated 

5.  imer&ce  elements  to  be  left  inactive  but  to  be  placed  between  elements  of  like 
nutterials. 

Column  Variable _ Explanation _ 

1-S  LEL(N,1)  Element  number  of  first  dement  to  be  placed 

6-10  LEL(N,2)  Material  number  of  the  element 

Information  for  the  next  7  placed  dements  is  supplied  in  the  next  14  five-column  fields. 

(3)  THIRD  CARD  (NODE  NUMBERS)  FORMAT(16IS) 

Nodal  point  numbers  of  16  nodes  to  be  asdgned  zero  displacement  can  be  supplied  on 
one  card.  A  maximum  of  SO  can  be  specified  in  one  loading  step.  Nodal  point  numbers,  N 
==  1  to  NONP,  are  to  be  specified  in  sequential  order. 
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Col«m0  Variable _ tkwi _ 

1-5  NP(N)  Nodal  point  mind)er  of  firtt  node  assigned  zeio  diq^laceinent. 


Infonnation  for  the  next  15  nodal  point  is  supplied  in  the  next  15  five-column  fi^ds. 

(c)  SEEPAGE 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  OR  KCS(N,3)  =  3.  Input  is 
handled  by  subroutine  SEEP.  Seepage  loads  are  determined  from  chides  in  pore  pressure 
^)ecified  as  DP(N)  on  nodal  point  cards,  or  fit>m  the  q)edfied  phreatic  levd  changes. 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT(15) 

Cdumn  Variable _ Explanation _ 

1-5  NCODE  Option  code  spediying  how  seqpage  is  input- 

s  0  If  specified  as  DP(N)  on  no^  point  cards  (section  9). 

-  1  If  to  be  calculated  using  the  new  phreatic  suifiice  input  on 
the  follovnng  cards. 

(2)  SECOND  CARD  (NUMBER  OF  PHREATIC  SEGMENTS)  FORMAT(IS) 
This  card  is  required  only  if  NCODE  =  1. 

Column  Variable _ Explanation _ 

1-5  NWAT  Number  of  phreatic  sur&ce  segment  and  points  used  to  specify 

the  new  phi^c  sur&ce.  NWAT  must  be  greater  than  or  equal 
to  2.  The  number  of  phreatic  sur&ce  segments  is  equal  to 
NWAT.  The  minimum  value  of  NWAT  is  30. 

(3)  THIRD  CARD  (PHREATIC  LEVEL  DATA)  FORMAT(6D10.2) 

Tins  card  is  required  only  ifNCODE  »  1. 

The  end  points  of  the  phreatic  surfiice  segments  delineating  the  new  and  old  phreatic 
surfiices  are  specified  as  x-coordinates  and  must  be  the  same  as  the  x-coordinates  of  a 
nodal  point  (Figure  A8).  Both  the  present  and  new  phreatic  surfiu^  are  assumed  to  be 
linear  between  the  boundary  x-coordiiutes.  The  left  hand  side  of  the  mesh  is  always  the 
first  boundary  x-coordinate  specified.  A  boundary  x-coordinate  on  the  old  phreatic  surfiice 
require,  usually,  qrecification  of  the  same  x-coordiru0e  on  the  new  phr^c  surfiice. 

Two  end  points  (x-coordinates),  with  assodated  new  and  old  phreatic  levds  (y- 
coordinates),  are  supplied  a  card  each.  All  end  points,  N  =  1  to  NWAT,  must  be 
specified. 
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Coiliaui  Variable 

I- 10  XW(N) 

II- 20  PREL(N) 

21-30  FUEL(N) 


_ Eiptonation  _ 

X-coordinate  boundaiy  the  levds  PREL(N)  and  FUEL(N)  on 
the  right  hand  side.  (Must  be  the  same  as  the  x-coordinate  of  a 
nodal  point.) 

Present  level  (y-coordinate)  of  the  phreatic  surface  at  XW(N)- 
New  level  (y-coordinate)  of  the  phreatic  surftce  at  XW(N). 


Information  for  the  next  end  point  is  supplied  in  the  next  3  ten-column  fields. 


Y 


Figure  A8.  Example  illustrating  phreatic  level  data  required 
for  the  seepage  loading/construction  case  (Clough  1984) 

(d)  DELETION  OF  BAR  ELEMENTS 

These  cards  supplied  only  if  KCS(N,1),  KCS(,N,2),  OR  KCS(N,3)  =4.  Input  is  handled 
by  the  nuin  program. 

The  deleted  bar  elements  remain  in  the  meidi  but  ^nth  zao  stififiiess.  The  force  the  bar 
element  carried  is  ^plied  to  the  fi'ee  node  or  nodes  at  its  ends. 
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This  loadiqg/construction  mode  cannot  be  specified  in  the  same  loading  as  fill  or 
concrete  placement. 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT(I5) 

Column  Variable _  Explanation 

1-5  NCARDS  Number  of  deleted  bar  elements.  There  is  no  limit  other  than 

the  number  of  bar  elements  presently  in  the  mesh. 

(2)  SECOND  CARD  (ELEMENT  NUMBERS)  FORMAT(16IS) 

Column  Variable  Explanation _ 

1*5  N  Element  number  of  bar  elements  to  be  deleted. 

Element  numbers  for  the  next  15  elements  are  supplied  in  the  next  15  five-column  fields. 

(e)  ADDITION  OF  BAR  ELEMENTS 

These  cards  are  suppUed  only  if  KCS(N,1),  KCS(N,2),  KCS(B,3)  =  5.  Input  is  handled 
by  the  main  program. 

Information  on  the  second  card  is  the  same  as  that  explained  for  the  Bar  Element  Card 
(section  12).  The  added  bar  elements  are  numbered  sequentially  form  NUMBAR  +  1, 
where  NUMBAR  is  the  number  of  bar  elements  in  the  mesh  before  the  present  loading 
step. 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT(I5) 

Coiunm  Variable _ Eipia»ation  _ 

1-5  NCARDS  Number  of  bar  elements  to  be  added.  Any  number  can  be 

added  in  a  loading  step,  however,  the  maximum  number  of  bar 
elements  that  can  be  in  the  mesh  (including  inactive  or  deleted 
elements)  is  15. 

(2)  SECOND  CARD  (ADDED  BAR  ELEMENT  DATA) 

FORMAT(4I5^D10.74D10.14>10.5) 


Column 

Variable 

Explanation 

1-5 

N 

Element  number  of  added  bar  element. 

6-10 

IB(N,1) 

Number  of  nodal  point  J 

11-15 

IB(N.2) 

Number  of  nodal  point  I 
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Column 

Variable 

Exidanation 

16-20 

IB(N,3) 

Spring  re^nse  type  code. 

21-30 

BAR(N,1) 

cos  a 

31-40 

BAR(N,2) 

ana 

41-50 

BAR(N,3) 

Prestress  force  in  the  bar  element. 

51-60 

BAR(N,4) 

Stiffiiess  of  bar  elemoit 

61-70 

BAR(N,5) 

Di^lacement  of  bar  element  at  activation 

Note  that  in  a  restart  analysis  the  added  bar  dement(s)  information  must  be  included  in  the 
bar  element  connectivity  cards  (Section  12)  and  the  value  for  NUMBAR  in  the  data 
control  card  (Section  2a)  must  be  updated. 

(0  BOUNDARY  PRESSURE  LOADING 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  OR  KCS(N,3)  =  6.  Input  is 
handled  by  subroutine  SURFLD.  Linear  pressure  distribution  are  assumed,  based  on  the 
pressure  specified  for  the  nodal  points. 


Column 

Variable 

Explanation 

1-5 

NLDS 

Number  of  loaded  boundaries.  There  is  no  limit  to  the  number 
of  loaded  boundaries  that  can  be  specified. 

6-10 

NOLDSX 

Number  of  loaded  boundaries  for  which  the  horizontal 
components  of  the  nodal  points  load  vectors  are  set  equal  to 

zero. 

11-15 

NOLDSY 

Number  of  loaded  boundary  for  wdiich  the  vertical  components 
of  the  nodal  points  load  vectors  are  set  equal  to  zero. 

(2)  SECOND  CARD  (LOADED  BOUNDARY  DATA)  FORMAT(2I5,4D10.2) 

Information  for  one  loaded  boundary  is  specified  on  each  card.  Nodes  I  and  J  are 
specified  counterclockwise  on  an  element  (Figure  A9).  A  normal  compressive  traction 
(pressure)  is  positive.  A  tangential  traction  (shear  stress)  is  positive  when  directed 
clockwise  (n^e  J  to  node  I)  as  shown  in  Figure  A9.  A  total  of  N  =  1  to 
(NLDS+NOLDSX+NOLDSY)  loaded  boundaries  must  be  spedfied.  NLDS  cards  are 
specified  first,  followed  by  NOLDSX  cards  and  the  NOLDSY  cards. 
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Column 

Variable 

Exotauation 

1-5 

I 

Nodal  point  luimber  of  the  first  imde  of  the  loaded  boundary. 

6-10 

J 

Nodal  point  numba  of  the  second  node  of  the  loaded 
boundary. 

11-20 

WSl 

Value  ofthe  normal  pressure  acting  at  node  I. 

21-30 

WS2 

Value  of  the  normal  pressure  actirig  at  node  J. 

31-40 

WS3 

Value  of  the  shear  stress  at  node  I. 

41-50 

WS4 

Value  of  the  shear  stress  at  node  J 

(g)  TEMPERATURE  LOADING 

No  cards  are  required  for  this  loading/constniction  mode.  Temperature  Card  (section 
7).  If  KCS(N,1),  KCS(,N,2),  or  KCS(N,3)  =  7,  then  the  values  of  DP(N)  are 
acknowledged  by  the  main  program  and  processed  as  ten4)erature  changes.  Note  that 
DP(N)  can  also  be  used  to  input  phreatic  level  changes  for  the  seepage 
loading/construction  mode.  Thus,  if  seepage  is  specified  as  being  input  through  values  of 
DP(N),  seepage  and  temperature  loading  cannot  be  included  in  the  same  loading  step. 
GenenJly  temperature  loading  requires  a  restart  analysis,  with  the  DP(N)  values  being 
changed  to  reflea  the  temperature  changes  prior  to  the  analysis. 

The  temperature  scale  used  (‘*C  or  **?)  must  correspond  to  the  coefiBcient  of  thermal 
expansion  designated  on  the  Material  Property  Cards  (section  5).  Temperature  changes 
are  typically  designated  for  structural  materials  only. 

(h)  CONCENTRATED  FORCE  OR  DISPLACEMENT  LOADING 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  =  8  (nodal 
displacements)  or  9  (nodal  loads).  Input  'is  handled  by  the  main  program.  Refer  to  the 
Loading  Information  Card  (Seaion  S)  for  instructions  regarding  loading/construction 
modes  8  and  9. 

(1)  FIRST  CARD  (CONTROL  DATA)  FORMAT(IS> 

Column  Variable _ Explanation _ 

1-S  NUMNDE  Number  of  loaded  or  displaced  nodes.  There  is  no  limit  to  the 

number  of  loaded  or  displaced  nodes  that  can  be  spedfied. 
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Figure  A9.  Exanq}le  of  irgmt  parameters 
for  positive  boundary  pressure  loading 
(Ebeling,  Peters  and  Chugh,  J992) 


3 


(2)  SECOND  CARD  DISPLACEMENT  DATA  FORMAT(2(15^(lX,lD14.7))) 
(LOAD  DATA]  FORMAT(2(I10,2D10.2)) 


Information  for  two  loaded  or  displaced  nodal  points  is  supplied  on  each  card.  A  total 
of  N  =  1  to  NUMNDE  nod^  must  be  specified.  Sign  convention  is  positive  to  the  right 
(potitive  x-direction)and  positive  up  (potitive  y-direction).  Nodal  points  specified  as  being 
loaded  or  displaced  do  not  have  to  be  in  numerical  order. 


Column  Variable 


1-5 

I 

7-20 

XI 

22-35 

Y1 

36-40 

J 

42-55 

X2 

57-70 

Y2 

_ Explanation _ 

Node  number  of  the  first  displaced  or  loaded  node. 

Component  of  displacement  or  force  in  the  x-direction  at  node 
I. 

Component  of  displacement  or  force  in  the  y-direction  at  node 

I. 

Node  number  of  the  second  displaced  or  loaded  node. 

Component  of  displa^ment  or  force  in  the  x-direction  at  node 

J. 

Component  of  displaconent  or  force  in  the  y-direction  at  node 
J 


If  there  is  no  second,  or  J  node  to  be  specified  on  the  last  card,  then  leave  the  second 
set  of  columns  blank. 


Displaced  nodes  are  to  be  included  in  Boundary  Condition  Codes  (Section  10) 

(i)  ELEMENT  MATERIAL  TYPE  CHANGE 

These  cards  are  supplied  only  if  KCS(N,1),  KCS(N,2),  or  KCS(N,3)  =10.  Input  is 
handled  by  the  main  program. 

The  material  type  of  the  specified  element  is  changed  before  the  analysis  of  the  loading 
step  which  specifies  the  change.  The  material  type  change  includes  modifying  the  values  of 
modulus,  E,  and  Poison's  ratio,  GUI,  and  zeroing  the  stresses,  SIG^,1).  Thus,  if  a 
material  type  change  is  specified  in  conjunction  with  boundary  loading,  in  the  same  loading 
step,  the  elements  whose  material  type  is  changed  will  respond  to  the  loading  with  new 
nuterial  properties. 

As  included,  this  loading/construction  mode  is  intended  to  physically  represent  the 
grouting  of  an  anchor.  At  a  given  step  in  the  analysis,  the  material  types  of  soil  elements 
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can  be  changed  to  rqMnesent  the  assumed  linear  elastic  grout  zone.  If  there  is  a  need  to 
change  material  types  for  ai^  other  reastm,  tins  can  be  dtme  by  stopinog  the  analysis  after 
the  apf»t^»iate  load  stq),  modifying  the  material  types  on  the  Matoial  Type  Doignation 
Card  (section  14D),  then  restarting  the  anafy^. 

(1)  FIRST  CARD  (CONTROL  CARD)  FORMAT(I5) 

The  maximum  number  of  elements  v^se  material  type  number  can  be  changed  in  a 
load  step  is  120.  The  excavation  and  material  type  change  loading/construction  modes 
cannot  be  specified  in  the  same  loading  step  since  the  same  variable,  LUL(N,I),  is  used  to 
input  data  for  both. 

Column  Variable _ Explamition _ 

1-5  NELCH  Number  of  donents  udmse  material  type  numba-  is  being 

changed. 

(2)  SECOND  CARD  (ELEMENT  DATA)  FORMAT(16IS) 

The  element  numbos  and  tww  material  type  numbers  of  8  dements  can  be  supplied  on 
one  card.  A  total  of  N~1  to  NELCH  dements  must  be  specified. 

Column  Variable _ EurtanatioB _ _ 

1-5  LUL(N,  1)  Element  number  of  first  dement  with  a  spedfied  new  material 

type  number. 

6-10  LUL(N,2)  New  nwterial  type  number  of  the  specified  dement. 

Information  for  the  next  seven  demoits  is  supplied  in  the  next  7  pairs  of  five-column 
fields. 
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APPENDIX  B:  SEQUENCE  OF  OPERATION 


1.  The  program  SOILSTRUCT  uses  the  direct  stiffiiess  method, 

solve  for  incremental  nodal  di^laconents  {ti}  resulting  from  iiKS’emental  loads  ^)plied 
to  the  nodal  points  {F}.  The  local  element  stiffiiess  matrices  are  first  formulated,  then 
assembled  into  the  global  stiffiiess  matrix.  Equivalent  nodal  loads  due  to  construction 
or  applied  loadings  are  assembled  in  the  incremental  load  vector,  {F}.  The  computed 
incremental  displacemems  are  then  used  to  compute  the  incremental  change  in  stress 
acting  at  the  center  of  the  elements.  The  values  of  total  stress  are  updated  by  the 
computed  incremental  changes  in  stress.  Tlw  total  stresses  are  then  used  to  revise  the 
elastic  moduli  used  in  the  formulation  of  the  element  stiffiiesses.  These  procedures  are 
repeated  for  each  iteration,  and  in  turn  for  each  substep  of  each  load  case  during  the 
analysis.  The  determination  of  initial  gravity  stresses  is  accomplished  in  a  similar 
fashion  and  can  be  viewed  as  an  initial  load  step,  with  the  nodal  loads  equal  to  the 
body  forces. 

2.  The  following  is  a  listing  of  the  names  of  each  of  the  33  subroutines  comprising 
the  program  SOILSTRUCT  and  a  brief  description  of  their  purpose: 

•  Main  Program.  The  main  program  serves  to  control  the  execution  of 
SOILSTRUCT.  It  calls  subroutines,  prints  input  data,  load  case  information, 
material  properties,  node  and  element  data,  aiul  boundary  conditions.  The  input 
data  for  excavation,  seepage,  embankment  construction,  and  boundary 
loadings  are  printed  in  their  respective  subroutines.  Calculated  equivalent  nodal 
loads  due  to  installation  or  deletion  of  bar  elements,  and  concentrated  forces  or 
displacements,  are  added  in  the  main  program. 

•  DETNA.  Subroutine  DETNA  calculates  the  number  of  degrees  of  freedom, 
determines  the  location  of  the  diagonal  terms  of  the  global  stiffiiess  matrix  in 
the  vector  SN,  and  computes  the  required  size  of  SN. 

•  INTTAL.  Subroutine  INTTAL  calcula  s^  ^^nd  prints  initial  stresses  for  a  gravity 
tum-on  analysis.  This  is  done  by  se(.r  ‘enr  al  calls  to  STRSTF,  OPTSOL,  and 
STRESS.  If  a  restart  analysis  is  ^ledfied,  INTTAL  reads  the  continuation  data 
and  initializes  the  material  property,  stress,  and  displacement  arrays  for  the 
nodes  and  elements. 

•  SUBSTP.  Subroutines  SUBSTP  controls  the  analytis  of  each  load  case  when 
substeps  are  specified  for  that  load  case.  SUBSTP  divides  the  calculated 
equivalent  nodal  point  loads,  applied  displacements,  or  temperature  changes 
imo  the  specified  number  of  equal  incronents  prior  to  performing  the  analysis. 

•  STRSTF.  The  terms  of  the  global  stiffiiess  matrix  are  assembled  in  subroutine 
STRSTF  by  sequential  calls  to  QUAD,  JTSTF,  BAREL  and  BEMSTF. 
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OPTSOL.  Subroutine  OPTSOL  solves  the  series  of  simultaneous  equations 
using  Crout  reduction  to  obtain  the  incremental  displacements. 

STRESS.  Subroutine  STRESS  conqnjte  stresses  and  strains  for  the  two- 
dimensional  elements  and  print  results.  STRESS  caU  MODCAL,  BAREL,  and 
JSTRES,  used  to  update  tlw  modulus  values  for  two-dinwnsional  elemoits, 
inter&ce  stiffiiesses,  and  bar  stiffiiesses  for  use  in  the  next  itovtion  or  load 
case. 

QUAD.  Subroutine  QUAD  compute  the  local  dement  stififiiess  matrix  and 
stress-strain  matrix  for  two-dimensional  elements,  computes  equivalent  nodal 
point  forces  due  to  temperature  changes  of  non-soil  dements,  and  compute 
equivalent  nodal  point  forces  due  to  gravity  forces. 

BAREL.  Subroutine  BAREL  compute  the  stiffiiess  of  the  bar  dements,  and 
updates  the  bar  forces. 

JTSTF.  Subroutine  JTSTF  computes  the  stifBness  of  the  inter&ce  elements. 

JSTTRES.  Subroutine  JSTRES  computes  the  stresses  and  relative 
displacements  for  the  inter&ce  elements,  updates  the  inter&ce  stiffiiess  values, 
and  prints  results. 

MODCAL.  Subroutine  MODCAL  updates  the  modulus  values  assigned  to  the 
soil  elements. 

BUILD.  Subroutine  BUILD  computed  the  nodal  point  loads  which  are 
equivalent  to  the  wdght  of  the  elements  representing  a  newly  placed 
embankment  lift  and  establishes  the  initial  stresses  and  material  properties  for 
these  newly  placed  elements. 

EXCAV.  Subroutine  EXCAV  calculates  the  stresses  acting  on  an  excavation 
boundary. 

EQNDFO.  Subroutine  EQNDFO  converts  the  stresses  calculated  by  EXCAV 
to  equivalent  nodal  point  forces,  i^di  are  added  to  the  incremental  load 
vector  by  EXCAV. 

SURFLD.  Subroutine  SURFLD  computes  equivalent  nodal  point  forces  due 
to  a  boundary  pressure  loading  iq)plied  along  the  face  of  an  elernem  and  adds 
these  computed  forces  to  the  incremental  load  vector. 

SEEP.  Subroutine  SEEP  calculates  equivalem  nodal  point  forces  due  to 
changes  in  the  phreatic  sur&ce  and  adds  these  forces  to  the  incremental  load 


vector.  The  nodal  p(^  forces  are  formulated  baaed  upon  changes  in  pore 
water  pressures. 

AUXOUT.  SulHXMitine  AUXOUT  writes  continuatimi  data  to  a  file  for  use  in 
subsequent  analyses. 

PRNCIP.  Subroutine  PRNCIP  calculates  principal  stresses  and  the  maximum 
shear  strain  for  two-dimensional  elements. 

PRNTFD.  All  non-zero  values  of  the  incremental  load  vector  are  printed  by 
subroutine  PRNTFD. 

GETFIL.  Subroutine  GETFIL  initializes  the  execution  of  SOILSTRUCT  by 
requesting  the  names  of  the  input  and  output  files  and  the  corresponding 
opening  and  closing  of  the  disc  storage  devices. 

NOTENS.  Subroutine  NOTENS  checks  for  tension  fiulure  within  interfoce 
elements,  reduces  excessive  stresses  within  fiuled  interfiice  elemoits  and 
removes  fiuled  elements  by  assigning  zero  stiffiiess. 

UNBALS.  Subroutine  UNBALS  conqxites  the  unbalanced  force  for  fiuled 
elonent  and  places  these  forces  in  the  n>  array  as  well  as  updates  the  dement 
stresses. 

PRNTFO.  Subroutine  PKNTFO  (Hints  the  contents  of  the  FO  array. 

PRNTJT.  Subroutine  PRNTJT  prints  the  contents  of  the  FO  array  for  all 
inteifoce  elements. 

NEWMOD.  Subroutine  NEWMOD  asrigns  the  new  moduli  to  all  dements. 

INTERF.  Subroutine  INTERF  diedcs  interfoce  dements  for  |)ossible  tension 
fiulure  at  the  center  of  the  dement  and  conqxites  the  fraction  of  iq>plied  forces 
resulting  in  zero  overshoot  normd  fcxce  for  that  interfiu:e  dement.  With  at  - 
0,  iq)plication  of  the  alpha  method  results  in  zero  noimal  stress  at  the  center 
line  of  the  imerfoce  dement. 

ALFA2D.  Subroutine  ALFA2D  locates  the  2-D  soil  dements  whose  stress 
levds  are  greater  than  1.00  and  conqxites  the  fraction  of  the  increment  in 
stresses  that  result  in  a  stress  levd  just  equal  to  1.00. 

COLAPS.  Subroutine  COLAPS  chedcs  to  if  excesrive  movement  has 
occurred. 


BEMSTF.  Subroutine  BEMSTF  calculates  the  equivalent  stiffiiess  matrix  of 
the  boundary  dement  system. 

MATRX.  Subroutine  MATRX  calculates  the  matrix  of  influence  coefficient  of 
the  boundary  donent  system. 

FUNC.  Subroutine  FUNC  computes  the  boundary  integral  necessary  in  the 
computation  of  the  matrix  of  influence  coeffidents. 

MELAN.  Subroutine  MELAN  computes  the  stresses  and  displacements  in  the 
Melan  fundamental  solution. 

SETCON.  Subroutine  SETCON  computes  the  constants  used  in  the 
computation  of  the  Melan  fundamental  solution. 

NMATRX.  Subroutine  NMATRX  computes  the  irtterpolation  matrix  N  for 
converting  nodal  tractions  to  nodal  forces. 

STFSYM.  Subroutine  STFSYM  extracts  the  symmetric  part  of  the  boundary 
element  stiffiiess  nuitrix. 
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User's  Guide:  Computer  Program  wHh  Interactive  Graphics  for 
Analysis  of  Plane  Frame  Structtses  (CFRAME) 

Survey  of  Bridge-Oriented  Design  Software 

Evaluation  of  Con^xjter  Programs  for  the  Design/Analysis  of 
Highway  and  Railway  Bridges 

User's  Guide:  Compiler  Program  for  Desi^i/Review  of  Curvi¬ 
linear  Conduits/Culverts  (CURCON) 

A  Three-Dimensional  Finite  Element  Data  Edit  Program 

A  Three-Dimensional  Stabflity  Anaiysia/Design  Program  (3DSAD) 
Reporti:  General  Geometry  Module 
Report3:  General  Analysis  Module  (CGAM) 

RefKxt  4:  Spedal-Purpoae  Modules  for  Dams  (CDAMS) 

Basic  User’s  Guide:  Computer  Program  for  Design  and  Analysis 
of  Inverted-T  Retaining  WaHs  and  FtoodwaHs  (TWDA) 

User's  Reference  Manual:  Computar  Program  for  Design  and 
Analysis  of  invertod-T  Retaining  Walls  and  FtoodwaHs  (TWDA) 

Documentation  of  Rnito  Element  Analyses 
Reporti:  Longview  Outlet  Works  Conduit 
Re^  2:  Anchored  Wall  Monolilh,  Bay  Springs  Lock 

Basic  Pile  Group  Behavior 

User’s  Guide:  Computer  Pro(pam  for  Design  and  Analysis  of  Sheet 
Pile  WaHs  by  Classicai  Methods  (CSHTWAI4 
Reporti:  Computattonal Processes 
Report  2:  Interactive  Graphics  Options 

Validation  Report  Computer  Proi^am  for  Design  and  Analysis  of 
Invertod-T  Retaining  W^  and  FtoodwaHs  (TWDA) 

User's  Guide:  Computer  Program  for  Design  and  Analysis  of 
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User's  Guide:  Computer  Program  for  Design  or  Invesdgation  of 
Orthogonal  Culverts  (CORTCUL) 

User'sGuide:  Computar  Program  for  Three-Dimensional  Analysis 
of  Building  Systems  (CTABS80) 

Theoretical  Basis  for  CTABS80:  A  Computer  Program  for 
Three-Dimensional  Analysis  of  BuBdlng  Systems 

User’sGuide:  Computar  Pro(^am  for  Analysis  of  Beam-Column 
Sbuctures  with  Nonlinear  Supports  (CBEAMC) 
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of  Shaloar  Foundations  (CKAR) 

User's  Quids:  Computer  Program  with  Interactive  Graphics  for 
Analysis  of  Plane  Frame  Structures  (CFRAME) 

User'sQuide:  Computer  Program  for  Generation  of  Engbteering 
Geometry  (SKETCH) 

User's  Guide:  Computer  Prof^am  to  Calculate  Shear,  Moment, 
and  Thrust  (CSMT)  from  Stress  Reeults  of  a  Two-Dimensional 
Finite  Elem^  Analysis 

Basic  Pile  Group  Behavior 

Rafeienoe  Manual:  Computy  Graphics  Program  for  Generation  of 
Engineering  Geometry  (SKETCH) 

Case  Study  of  Six  Major  General-Purpose  Finite  Element  Programs 

User'sQuide:  Computer  ProfFam  for  Optimum  Dynamic  Design 
of  Nonlinear  Metal  Plates  Under  Blast  Loading  (CSDOOR) 

User's  Guide:  Computer  Program  for  Detennining  Induced 
Stresses  and  Consolidation  Settiements  (CSETT) 

Seepage  Analysis  of  Confined  Flow  Problems  by  the  Method  of 
Fragments  (CFRAG) 

User's  Guide  for  Computer  Program  CGFAG,  Concrete  General 
Flexure  Analysis  with  Graphics 

Computer-Aided  Drafting  and  Design  for  Corps  Structural 
Engineers 

Dedsion  LogicTsbIe  Fomudaflon  of  ACI 318-77,  Building  Code 
Requirements  for  Reinforeed  Concrete  for  Automated  Con¬ 
straint  Processing.  Volumes  I  and  II 

A  Case  Committee  Study  of  Finite  Element  Analysis  of  Concrete 
Flat  Slabs 

User’sQuide:  Computer Prosyam for Two-Dimensionai Analysis 
of  U-Frame  Sfructures  (CUFRAM) 

User'sQuide:  For  Concrete  Strentyi  investigation  and  Design 
(CASTR)  in  Accordance  with  ACI 318^ 

RnRe-Element  Method  Package  for  Sdvfrig  Steady-State  Seepage 
Problems 

User's  Guide:  A  Three  Dbnensional  StaMIty  Analysis/Design 
Program  (3DSAO)  Module 
Report  1:  Revision  1:  General  Geometry 
Reixirt2:  General  Loads  Module 
Reports:  Free  Body  Moctole 
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Finite  Element  Studtes  of  a  Horizontally  Framed  Miter  Gate 

Report  1:  Initial  and  Refined  Finite  Element  Models  (Phases 
A,  B.  and  C),  Volumes  I  and  II 
Report  2:  Sim^ified  Frame  Modal  (Phase  D) 

Reports:  Altsmale Configuration MHer Gals Finils Element 
Sturfies-Open  Section 

Report4:  AUemate  Configuration  MHer  Gale  Finils  Element 
Studies-dosed  Sections 

Reports:  Aitsmate Configuration  Mitar Gale  Rnto  Element 
Studtos-Adcfltional  Closed  Sections 
Reports:  Elastic  Bucking  of  Girders  in  Horizontally  Framed 
Mhsr  Galas 

Report  7:  A, -  .iicationandSuniiTtary 


User's  Guide:  UTEXAS2  Slope-Stabilty  Pack^je;  Volume  I. 
User's  Manual 


SIkfing  Stability  of  Concrete  Structures  (CSUDE) 

Critaria  Specifications  for  and  Vaidation  of  a  Computer  Program 
for  the  Design  or  Investigation  of  Horizontally  Fiained  Miter 
Gales  (CMITER) 

Procedure  for  State  Analysis  of  Gravity  Dams  Using  the  Rnite 
Element  Method  -  Phase  la 

Usor’sGuMe:  Computer  Program  for  Analysis  of  Planar  Grid 
Structures  (CGRID) 

Development  of  Design  Formulas  for  RBibed  Mat  Foundations 
on  Expansive  Sois 

User's  Guide:  Pile  Group  Graphics  Display  (CPGG)  Post- 
prccessorto  CPGA  Progtwn 

User’s  Guide  for  Design  and  investigation  of  Horizontally  Framed 
MHer  Gates  (CMITER) 

User's  Guide  for  Revised  Computer  Program  to  Calculate  Shear, 
Moment  and  Thrust  (CSMT) 

User’s  Guide:  UTEXAS2  Slope-Stabiiity  Package;  Volume  II. 
Theory 

User’s  Guide:  PHe  Group  Analysis  (CPGA)  Computer  Group 

CBASIN-Structural  Desigft  of  Saint  Anthony  Fans  Staling  Basins 
Aooonfing  to  Corps  of  Engineers  Criteria  for  Hydraulc 
Structures;  Computer  Proi^am  X0098 
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CCHAN-Structurai  Design  of  Rectangular  Channels  Accorrfing 
to  Corpe  of  Engineers  Criteria  for  Hydraulic 
Structures;  Computer  Program  X0097 

The  Response-Spectrum  Dynamic  Analysis  of  Gravity  Dams  Using 
the  Rnite  Element  Method;  Phase  II 

State  of  the  Art  on  Expert  Systems  Applications  in  Design, 
Construction,  and  Maintenance  of  Structures 

User’s  Guide:  Computer  Program  fw  Design  and  Analysis 
of  Sheet  Pile  Walls  by  Classical  Methods  (CWALSHT) 

Investigation  and  Design  of  U-Frame  Structures  Using 
Program  CUFRBC 

Volume  A:  Program  Criteria  and  Documentation 
Volume  B:  User’s  Guide  for  Basins 
Volume  C:  User’s  Guide  for  Channels 

User’s  Guide:  Computer  Program  for  Two-Disnensional  Analysis 
of  U-Frame  or  W-Frame  Structures  (CWFRAM) 

User’s  Guide:  Pile  Group-Concrete  Pile  Analysis  Program 
(CPQC)  Preprocessor  to  CPGA  Program 

Application  of  Finite  Element,  Grid  Generation,  and  Scientific 
Visualization  Techniques  to  2-D  and  3-D  Seepage  and 
Groimdwater  Modeling 

User’s  Guide:  Computer  Program  for  Design  and  Analysis 
of  Sheet-Pile  Walls  by  Classical  Methods  (CWALSHT) 

Including  Rowe's  Moment  Reduction 

User’s  Guide  for  Concrete  Strength  investigation  and  Design 
(CASTR)  in  Accordance  with  ACI 318-89 

Finite  Element  Modeling  of  Welded  Thick  Plates  for  Bonneville 
Navigation  Lock 

Introduction  to  the  Corr^xitation  of  Response  Spectrum  for 
Earthquake  Loading 

Concept  Design  Example,  Computer  Akteu  Structural 

Modeling  (CASM) 

Reporti:  SchemeA 
ReiMrt  2:  Scheme  B 
Reix)rt3:  Scheme  C 

User's  Guide:  Computer-Aided  Structural  ModeOng 
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Optimization  of  Steel  Pile  Foundations  Using  Optimality  Criteria  Jun  1992 
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Knowledge-BasedExpertSystemforSelectionandDesign  Sep  1992 
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Evaluation  of  Thermal  and  Incremental  Construction  Effects  Sep1992 

for  Monoliths  AL-3  and  AL-5  of  the  Melvin  Price  Locks 
and  Dam 
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